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Abstract. We study the thermodynamics of the Hamiltonian Mean Field (HMF) model with an external 
potential playing the role of a "magnetic field" . If we consider only fully stable states, this system does not 
present any phase transition. However, if we take into account metastable states (for a restricted class of 
perturbations), we find a very rich phenomenology. In particular, the system displays a region of negative 
specific heats in the microcanonical ensemble in which the temperature decreases as the energy increases. 
This leads to ensembles inequivalence and to zeroth order phase transitions similar to the "gravother- 
mal catastrophe" and to the "isothermal collapse" of self-gravitating systems. In the present case, they 
correspond to the reorganization of the system from an "anti-aligned" phase (magnetization pointing in 
the direction opposite to the magnetic field) to an "aligned" phase (magnetization pointing in the same 
direction as the magnetic field). We also find that the magnetic susceptibility can be negative in the mi- 
crocanonical ensemble so that the magnetization decreases as the magnetic field increases. The magnetic 
curves can take various shapes depending on the values of energy or temperature. We describe hysteretic 
cycles involving positive or negative susceptibilities. We also show that this model exhibits gaps in the 
magnetization at fixed energy, resulting in ergodicity breaking. 

PACS. 5.20.-y Classical statistical mechanics - 05. 45. -a Nonlinear dynamics and chaos - 05.20.Dd Kinetic 
theory - 64.60.De Statistical mechanics of model systems 
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1 Introduction 

The statistical mechanics of systems with long-range in- 
teractions has recently been the object of an intense ac- 
tivity |ll2l3l4j . In particular, the Hamiltonian Mean Field 
(HMF) model |5l6l7l8l9ll0in] . which is a prototype for 
systems with long-range interactions, has been particu- 
larly studied by statistical mechanicians in order to illus- 
trate the beauty and the richness of these systems. The 
HMF model can be viewed as a XY spin system with in- 
finite range interactions or as a one dimensional model of 
particles moving on a ring and interacting via a long-range 
potential truncated to one Fourier mode (cosine poten- 
tial) . In that second interpretation, it shares many analo- 
gies with self-gravitating systems |7l8l9ll0fTT| but is much 
simpler to study since it avoids difficulties linked with the 
singular nature of the gravitational potential at the origin 
and the absence of a natural confinement |12I13I14] . 

The HMF model has been shown to possess many 
interesting features both at equilibrium and out-of- 
equilibrium. The dynamical evolution of the HMF 
model in the microcanonical ensemble displays non 
Boltzmannian quasistationary states (QSS) |10I15I16] 
that are stable steady states of the Vlasov equation 
|7lllll7ll8ll9l20l21l22j . out-of-equilibrium phase transi- 
tions explained by Lynden-Bell's theory of violent re- 
laxation |23l24l25l26l27l28l29j or by dynamical processes 



[30131132155] . re-entrant phases |24I29| . negative kinetic 
specific heats |34I35) , incomplete violent relaxation |19I35) , 
vanishing of the Lyapunov exponent |36l37j . aging [38], 
glassy dynamics [32], collective oscillations [JO]: non- 
ergodicity [41142] . dynamical phase transitions [33], alge- 
braic velocity correlation functions and anomalous diffu- 
sion [44I45I46I47I48I49I50] . front structure of the velocity 
distribution tails |51| and slow convergence towards the 
Boltzmann distribution [15116118145152] explained by ki- 
netic theory [4153154] . The influence of an external noise 
and the coupling with a thermal bath (canonical ensem- 
ble) have also been considered [Tll5 5 56 5 7135] . as well 
as the effect of collisions between particles [SS]. At equi- 
librium, the caloric curve f3{E) displays a second order 
phase transition between a homogeneous phase at high 
energies/temperatures and a clustered phase at low en- 
ergies/temperatures [5I6I7I8I9I10I11I60I61] . This is similar 
to the collapse (or Jeans instability) of self-gravitating sys- 
tems [7'9'11'21J. However, contrary to self-gravitating sys- 
tems (12..13.14j , the caloric curve does not exhibit negative 
specific heats and the ensembles are equivalent [3]. 

New features arise when the system is submitted to an 
external potential playing the role of a magnetic field. The 
magnetic field imposes a specific direction to the magne- 
tization and breaks the rotational symmetry of the orig- 
inal HMF model. In this paper, we study in detail the 
thermodynamics of this system. If we consider only fully 
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stable states, there is no phase transition. However, if 
we take into account metastable states (for a restricted 
class of perturbations), we find a very rich phenomenol- 
ogy. In particular, the system displays a region of nega- 
tive specific heat in the microcanonical ensemble in which 
the temperature decreases as the energy increases. This 
leads to ensembles inequivalence and to zeroth order phase 
transitions similar to the "gravothermal catastrophe" and 
"isothermal collapse" of self-gravitating systems. In the 
present case, they correspond to the reorganization of 
the system from an "anti-aligned" phase (magnetization 
pointing in the direction opposite to the magnetic field) to 
an "aligned" phase (magnetization pointing in the same 
direction as the magnetic field) . We also find that the mag- 
netic susceptibility can be negative in the microcanonical 
ensemble so that the magnetization decreases as the mag- 
netic field increases. The magnetic curves can take various 
shapes depending on the values of energy or temperature. 
We describe hysteretic cycles involving both positive and 
negative magnetic susceptibilities. We also show that this 
model exhibits gaps in the magnetization at fixed energy, 
resulting in ergodicity breaking. 

A preliminary study of this model has been performed 
by Velazquez & Guzman |62j who derived several inter- 
esting results. We believe, however, that their study car- 
ries out some mistake^ and that complements are needed. 
This is the object of the present paper. Our paper closely 
follows the presentation and the methodology exposed in 
|lll61j for the usual HMF model without magnetic field. 



2 Thermodynamical approach 

2.1 The HMF model with a magnetic field 

The HMF model is a system of N particles of unit mass 
m = 1 moving on a circle and interacting via a cosine 
potential. When an external potential (named "magnetic 
field") H is imposed, the dynamics of these particles is 
governed by the Hamilton equations 



de, _ an 

dt dvj ' 



dvi 



on 



^ = 2 ^ ~ 4^ ^ - ) - i/^ cos0„ (1) 

i—l i^j i—1 

where 0i G [— tt, tt] and —oo < Vi < +oo denote the po- 
sition (angle) and the velocity of particle i and k is the 
coupling constant (we assume here that fc > 0). We have 
assumed, without loss of generality, that the field H is 
directed along the x-axis. Like for the usual HMF model, 
the proper thermodynamic limit corresponds to N ^ +oo 
in such a way that the rescaled energy e = SirE/kM'^ and 



^ We find, indeed, different conditions of stability. We think 
that the reason comes from the fact that Velazquez & Guzman 
|62j work in terms of the internal energy, while we work in 
terms of the total energy taking into account the contribution 
of the magnetic field. 



the rescaled inverse temperature 77 = pkM/Air remains 
of order unity [11]. We can take k ~ 1/A^ [5] which is 
the Kac prescription. In that case, the energy is exten- 
sive {E/N r-^ 1) and the temperature intensive (/3 ~ 1) 
but the system remains fundamentally non-additive [1]. 
For TV — !■ +00, the mean field approximation is exact 
[5l63j and the iV-body distribution function is a prod- 
uct of -V one-body distributions: P]s[{Oi^ Vi, ...^d]y,VN,t) = 
Pi{9i,vi,t)...Pi{9N,VN,t). 

Let us introduce the distribution function fiO^v^t) = 
NPi{6, V, t). For a fixed interval of time and N +00, the 
evolution of the distribution function f{9, v, t) is governed 
by the Vlasov equation [63] : 

at de ae av ae av ' ^ ' 



where 



2tt 



27r 



cos{9 - 9')p{e',t)d9', (3) 



is the self-consistent potential generated by the density of 
particles p{0,t) = J f{d,v,t)dv and 



(4) 



is the external potential. The mean force acting on a par- 
ticle located in 9 is {F){e,t) = -d<P/d0{e,t) - <P'^^t{e). 
Expanding the cosine function in equation ([3]), we obtain 



,t) = —Bx cos 9 — By sin 9, 



where 



Br 



B,, 



2-K 

k 

2^ 



p{9,t) COS 9 d9, 



p{9,t) sm9d9, 



(5) 



(6) 



(7) 



are proportional to the magnetization ^ J p cos 9 d9 
and by = jj J psm9d9. The magnetization can be viewed 
as the order parameter of the HMF model. 
Let us introduce the mass 



M[p] 



pd9, 



(8) 



and the mean field energy 

E[f] = ^ / ^^^^ + Ijp^d9 + J p^,,t d9, (9) 

where the first term is the kinetic energy, the second term 
the potential energy of interaction and the third term the 
potential energy due to the external field. Using equations 
(l5])-([7]), the potential energy can be expressed in terms of 
the magnetization as 



W 



-BrH. 



(10) 
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Wc also introduce the Boltzmann entropy 



S[f] 



/ln( ^ ) d9dv, 



and the Boltzmann free energy 

F[f] = E[f]^TS[f], 



(11) 



(12) 



where T = 1//3 > Ois the temperature. In the micro- 
canonical ensemble, the statistical equilibrium state is de- 
termined by the maximization problem 



max{S[f]\E[f] 



E, M[f] = M} . 



(13) 



In the canonical ensemble, the statistical equilibrium state 
is determined by the minimization problem 



mm I M[/] = M}, 



(14) 



The Boltzmann entropy functional (|lip and the maximum 
entropy principle (|13p can be justified by a standard com- 
binatorial analysis (see, e.g., [H]). The distribution func- 
tion f{9,v) that is solution of ([13]) is the most probable 
macroscopic state, i.e. the macrostate that is the most 
represented at the microscopic level, assuming that the 
accessible microstates (with the proper values of mass and 
energy) are equiprobable. 



2.2 Critical points 

We shall first determine the critical points of these vari- 
ational problems. The critical points of the maximization 
problem (|13p are determined by the variational principle 



SS - I35E - aSM = 0, 



(15) 



where /? = 1 /T and a are Lagrange multipliers associated 
with the conservation of energy and mass. The critical 
points of the minimization problem (jl4p are determined 
by the variational principle 



6F + aT6M = 0, 



(16) 



where a is a Lagrange multiplier associated with the con- 
servation of mass. Since T is fixed in the canonical ensem- 
ble, it is clear that equation (fT6|) is equivalent to equation 
Therefore, the optimization problems and (|14l) 
have the same critical points. Performing the variations in 
equations p3|) and , we find that the critical points are 
given by the mean field Maxwcll-Boltzmann distribution 



(17) 



where A' = Me"^"" is a constant and <l>tot{0) = ?^(6') + 
'^ext{0) is the total potential in 6. Integrating over the 
velocities, we get the mean field Boltzmann distribution 



p(0) =Ae-^*'-W, 



(18) 



where A — {2n/ PY^^A' . Using the expressions (H)) and 
^ of the potential, the distribution function (flT)) can be 
rewritten 



f{0,v) 



(19) 



In the following, we consider critical points whose mag- 
netization B is parallel to the magnetic field H so that 
By ~ and = B (we show in Appendix |B] that there 
is no critical point with By ^0). In that case, the forego- 
ing expression takes the form 



mt;)=A'e-4*-(^+^)-^t 
The corresponding density profile is 

p(0) =Ae^(^+^)™^^ 



(20) 



(21) 



The amplitude A and the magnetization B are determined 
by substituting equation (PT|) in equations and 
This yields 

A 



and 



27r/o(/3(B + i/))' 

2ttB ^ h{/3{B + H)) 
kM ~ Io{l3{B + H))' 



(22) 



(23) 



where In{x) is the modified Bessel function of order n. 
Equation (j23p determines the magnetization i? as a func- 
tion of the temperature T. Then, A is given by equation 
(221) ■ Finally, the distribution function and the density 
profile can be written 



1/2 



p{9)e 



p{9) = 



M 



27rIo{l3{B + H)) 



MB+H) cose 



(24) 



(25) 



where B is determined in terms of T by equation (j23p . 
The study of the self-consistency relation (j23[) will be per- 
formed graphically in Section] 



2.3 Thermodynamical parameters 

Let us now determine the expressions of the energy, en- 
tropy and free energy. For the Maxwcll-Boltzmann distri- 
bution (im, the kinetic energy is 



K = -MT. 
2 



(26) 



Combining this relation with equation (|10p . we find that 
the total energy E — K + W is given by 



1 ttB^ 27r 

E = -MT ; -BH. 

2 k k 



(27) 



The series of equilibria giving T as a function of E is 
determined by equations ([^5)) and ([?7)) by eliminating B. 
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Fig. 1. Series of equilibria (caloric curve) giving the inverse 
temperature r; as a function of the energy e for h = 0.01. In 
this figure, and in the following figures, we have indicated the 
values of x that parameterizes these curves. 
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Fig. 2. Magnetization b (order parameter) as a function of the 
inverse temperature 77 for h = 0.01. It exhibits a turning point 
of temperature at r]c{h). 



The relation that gives the magnetization _B as a function 
of the energy E is determined by equations ([^ and (P?]) 
by eliminating T. Finally, using equations ([TT|) and (IM|) . 
the entropy is given by 



5* = -MlnT- I plnpde, 



(28) 



up to a term + ^A/ ln(27r) + MlnAI. Using equation 
(l25l). it can be rewritten 



S 



iA/lnT + A/ln/o(/3(S + H)) - 



H) 



(29) 



up to a term + |A/ln(27r). The relation between the 
entropy S and the energy E is determined by equations 
(P7)) and (1121) by eliminating T and B. Using equa- 
tions ([^ and (P7)) . the free energy (fl^ is given by 



F = \mT - ImTIuT- MT\nIo[(i{B + H)) 



(30) 



up to a term -iil/T- |A/rin(27r). The relation between 
the free energy F and the temperature T is determined 
by equations (PU)) and ([25)) by eliminating B. 

It is convenient to write these equations in parametric 
form by introducing the parameter x = (3{B + H). Then, 
defining h = 2TrH/{kM), we obtain 



b = 



2ttB Ii{x) 



kM 

fikM 
47r 



&t:E 
kAP 



1 



2{b{x) + h)' 
- 2b{xf - Ahb{x) 



(31) 



(32) 



(33) 
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Fig. 3. Magnetization h (order parameter) as a function of the 
energy e for h = 0.01. It exhibits a turning point of energy at 
ec(/i). 
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Fig. 4. Entropy s as a function of the energy e for h = 0.01. 
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Fig. 5. Same as Figure |4] (zoom) 



r| (h) 




h = 0.01 



Fig. 6. Free energy / as a function of the inverse temperature 
77 for /i = 0.01. 
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Fig. 7. Same as Figure [6] (zoom) 



S 1 

s = — = —-lmj{x) +lii/o(a;) — xb{x), (34) 



M 



f 



eix) 



ri{x) 



s{x), 



(35) 



where unimportant constants have been omitted in the 
expression of the entropy. From these relations, we can 
obtain the curves T{E), B{T), B{E), S{E) and F{T) in 
parametric form. These curves are plotted in Figures [T][7] 



for h = 0.01. Wc have chosen a small value of h in order 
to show the connection with the results of the ordinary 
HMF model {h = 0) [11161] . Other curves are represented 
in Section [^21 The relations (PT|) - (|55|) characterize all the 
critical points of p3|) and (|14p . The structure of the series 
of equilibria will be studied in Section [2.4l Then, we must 
select among the critical points those that are (local) max- 
ima of S at fixed E and M (microcanonical ensemble) and 
those that are (local) minima of F at fixed M (canonical 
ensemble). This will be the object of Section [2761 



2.4 Series of equilibria 

By definition, the series of equilibria is formed by all the 
critical points of entropy at fixed mass and energy or, 
equivalently, by all the critical points of free energy at 
fixed mass (we have already indicated that the variational 
problems and (fH)) have the same critical points). The 
stability of these critical points, in each ensemble, will be 
investigated in Section [^751 bv using the Poincare theorem 
and in Sections [3] and |4] by studying the sign of the second 
order variations of entropy or free energy. In this section, 
we simply describe the structure of the series of equilibria. 

Let us assume that ft, > to fix the ideas (the case 
h < can be treated symmetrically). To understand the 
following discussion, it is required to look in parallel at 
Figures [T][7] and at the asymptotic expansions of Appendix 
Rl The graphical construction of Figure [T7] can also be 
useful. We need to distinguish two curves: (i) the series 
of equilibria with positive magnetization (same sign as 
the imposed field) corresponds to < a; < +00. It will 
be called the "aligned" phase. For x 0, we find that 
b ^ 0, 7] ^ and e — > +00. More precisely, i] ^ 1/e, 
b ^ hrj and b ^ h/e. For x — > +cxd, we find that 6 — > 1, 
•q — >■ +00 and e — > €min{h) = —2 — Ah. More pre- 
cisely, 77 ~ 2/(e - e„„„(ft)), 1 - 6 - 1/(4(1 + h)ri) and 
1 ~ b ^ (e — emin(^))/(8(l + h)); (ii) the series of equi- 
libria with negative magnetization (sign opposite to the 
imposed field) corresponds to —00 < x < xo{h). It will 
be called the "anti-aligned" phase. For x — !> —00, we find 
that b -1, rj +00 and e <^'min{^) = -2 + Ah. More 
precisely, ^ 2/(6 - e;„,„(ft)), -1 - 6 ^ -1/(4(1 - h)r,) 
and -1-6 ~ -(e-e;„„(/i))/(8(l-/i)). For x xo{h), we 
find that b 6o(/i) = —h, 7] — > +00 and e — > eo{h) = 2/i^. 
More precisely, 77 ~ l/(e — eo(/i)), bo{h) — b ^ ~xo{h)/{2ri) 
and bo{h)—b ^ —{xQ{h)/2){e—eo{h)). The series of equilib- 
ria with negative magnetization presents a turning point of 
temperature at 77 = r]c{h) corresponding to x = xj;""°(/i) 
and a turning point of energy at e = edh) corresponding 
to a; = a;™*'=™(/i). Note that a;™"°(/i) < a;™*'='^°(/i) so that 
the turning points {e{r]c{h)) , rjdh)) and (edh) , r]{ec{h))) 
differ. The "anti-aligned" phase exists only for ft, < 1. For 
h — > 1^, we find that 77c(ft) — > +00, ec{h) — > 2, eo{h) 2, 

CznW ^ 2 and 6o(/i) ^ -1- 

The asymptotic results described previously can be un- 
derstood easily. 

At very high temperatures T +00 (i.e. 77 — > 0) in the 
canonical ensemble or at very high energies e ^> +00 in 
the microcanonical ensemble, the long-range interaction is 
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negligible with respect to thermal motions and the system 
behaves essentially like a noninteracting perfect gas. The 
equilibrium state is a spatially homogeneous configuration 
(b ~ 0) with Maxwcllian velocity distribution. This is the 
unique (global) entropy maximum at fixed mass and en- 
ergy or the unique (global) free energy minimum at fixed 
mass. In that high energy /temperature limit, e ~ i/rj. 

Let us now consider the zero temperature state T = 
(i.e. r] +oo) in the canonical ensemble or the min- 
imum energy state in the microcanonical ensemble. To 
that purpose, we have to determine the minimum en- 
ergy at fixed mass. The kinetic energy is minimized by 
assigning the velocity t; = to each particle. Then, we 
have to minimize the potential energy w — —2{b'^ + 2bh). 
It is easy to show that the global energy minimum is 
f{0,v) = MS{v)6{6) corresponding to a magnetization 
b ~ +1 and an energy e,„i„(/i) = —2 — Ah. There also 
exists a local energy minimum f{6,v) = MS(v)d{0 — tt) 
corresponding to a magnetization b = —1 and an en- 
ergy i'„iinW — ~2 -|- Ah and a local energy maximum 
f{e,v) = ^M6{v)[S{e - Oq) 5{e + 0^)] with cos 6*0 = -h 
corresponding to a magnetization b = —h and an energy 
eo{h)=2h'^. 

For h ^ 0, we recover the results of [11I61[ valid for 
/i = 0. In particular, e„j„ = e^j„ = -2, cq = and Ec = 
rjc = 1. In that case, the points {^iilc) , r]c) and (ec,'7(ec)) 
coincide and the curve 77(e) forms a "spike" at {ec,r]c) = 
(1,1) leading to a second order phase transition [J]. On 
the other hand, the branch of inhomogeneous solutions 
becomes degenerate due to the invariance by rotation of 
the phase of the magnetization. 



2.5 Specific heat 

The specific heat is defined by C = dE/dT in both en- 
sembles. Using dimcnsionless variables, the specific heat 
per particle c = C/N can be written 



1 de 



2d{l/r,)- 

From equations ([5^ and ([55]) we easily obtain 
dil/v) 



dx 



K + -^'(^)' 

xri{x) X 



and 



dt 
dx 



1 



xrj{x) 

where b' ix) is given by 



2b' {x) 



1 



x rj{x) 



(36) 



(37) 



(38) 



(39) 



To obtain this expression, we have used the identities 
Iq{x) = Ii{x) and 
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Fig. 8. Specific heat c as a function of the inverse temperature 
rj for h = 0.01. We find that c — 1/2 for 77 — > (corresponding 
to X — >■ 0), c — )■ 1 for -q — +00 (corresponding to 2; -^■ ±00) 
and c — i> 1/2 for r/ — s- +00 (corresponding to a; xo{h)). In 
tire anti-aligned phase (—00 < x < xo{h)), the specific heat 
diverges like c oc ±(77 — ric{h))~^^^ — >■ ±co at the turning point 
of temperature 7?c(/i) (corresponding to a:: = x'^°'"°{h)). Close to 
ric(h), there is a region of negative specific heat (correspond- 
ing to xT'''°{h) < x < xT''"'''{h)). The study of Section \2B 
shows that the condition of instability in the canonical ensem- 
ble corresponds to x'^°''^°{h) < x < xo{h). In particular, the 
states with negative specific heats are unstable in the canoni- 
cal ensemble in agreement with general theorems of statistical 
mechanics. 



Therefore, the specific heat can be written 

_ 1 - 2b'{x) [t]{x) - x^] 
2 [1 - 2b'{x)T]{x)] 



(41) 



The specific heat is infinite at the turning point of temper- 
ature ric{h) corresponding to a;J;°"°(/i) and is zero at the 
turning point of energy edh) corresponding to x^^'^^°{h). 
It is represented as a function of the inverse temperature 
•q (appropriate to the canonical ensemble) in Figure [5] and 
as a function of the energy e (appropriate to the micro- 
canonical ensemble) in Figure [9] The series of equilibria 
with positive magnetization has positive specific heat for 
any < x < -t-00. Alternatively, there exists a region of 
negative specific heat in the series of equilibria with nega- 
tive magnetization between x'^^°{h) and x™*'^''°(/i) (while 
the specific heat is positive in the rest of the curve) . This 
region of negative specific heat can be seen directly on the 
series of equilibria 7y(e) of Figure [T] (see a zoom in Figure 
To]). In this region, the temperature decreases when the 
energy increases! This is the first occurrence of negative 
specific heats for the Boltzmannian statistical equilibrium 
state of the HMF modefl Indeed, in the absence of mag- 
netic field /i = 0, the specific heat is always positive but 
undergoes a discontinuity Ac = 2 at 7; = 77c or e = Ec jllj . 

It is a general result of statistical mechanics that stable 
states in the canonical ensemble have positive specific heat 
since the specific heat C ~ f3'^{{AE)'^) > measures the 



I'^ix) = In~lix) Inix). 

X 



(40) 



^ Out-of-equilibrium distributions 
heats have been found in [34I35| . 
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Fig. 9. Specific heat c as a function of energy e for h = 0.01. 
We find that c — > 1/2 for e — )■ +oo (corresponding to a; — >■ 0), 
c — >■ 1 for e — emin(^) and e ^> e^i„(/i) (corresponding to 
X- ^> ±cxd) and c — s> 1/2 for e — !> eo(ft) (corresponding to 
X — >■ 2;o(/i)). In the anti-aligned phase (— oo < a; < a:o(/i)), 
the specific heat diverges like c oc — (e — e(?7c(/i)))~^ — > icxD 
at the energy e{ric{h)) associated with the turning point of 
temperature r]c{h) (corresponding to a; = x'^°'"°{h)). Close to 
e(?7c(/i)), there is a region of negative specific heat (correspond- 
ing to xT"°(h) < X < xr'"'''{h)). The specific heat becomes 
zero at the turning point of energy €c{h) (corresponding to 
X = 2;J^""'°(/!,)). The study of Section [231 shows that the condi- 
tion of instability in the microcanonical ensemble corresponds 
to x^^°^°{h) < X < xo{h). Therefore, the states with negative 
specific heats are stable in the microcanonical ensemble. 



variance of the fluctuations of energy. Therefore, we can 
already conclude that the states between points CE and 
MCE in the series of equilibria (see Figure ITO)) are thermo- 
dynamically unstable in the canonical ensemble. However, 
negative specific heat is not a necessary condition of insta- 
bility in the canonical ensemble. On the other hand, stable 
states in the microcanonical ensemble can have negative 
specific heats [1]. Therefore, the study of the specific heat 
is not sufficient to settle the stability /instability of the 
system. 

2.6 Poincare theorem 

The stability of the critical points in each ensemble can be 
easily obtained by applying the Poincare theory of linear 
series of equilibria. This method has been extensively used 
in astrophysical problems (see, e.g., |13ll4j ). Linear series 
of equilibria are plotted in Figs. [T][7]for h = 0.01 and in 
Figs. [TTIfT3l for different values of h. 

Let us first consider the canonical ensemble. We have 
to determine the minima of free energy at fixed mass. A 
global minimum will be called fully stable (S), a local min- 
imum will be called metastable (M) and a maximum or 
a saddle point will be called unstable (U). To apply the 
Poincare theorem, we just have to plot e as a function of rj 
(see Figure [T] rotated by 90°). The series of equilibria with 
positive magnetization exists for any inverse temperature 
77 > 0. We know that the system is stable at high tempera- 
tures since it becomes equivalent to a classical gas without 
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Fig. 10. Enlargement of the series of equilibria Tj{e) in the re- 
gion of negative specific heats for h = 0.01. In the canonical 
ensemble, the system becomes unstable after the turning point 
of temperature (c = oo) and in the microcanonical ensemble, 
the system becomes unstable after the turning point of energy 
(c = 0). In particular, in the region of negative specific heats 
(i.e. between CE and MCE) the system is stable in the mi- 
crocanonical ensemble but unstable in the canonical ensemble. 
This can be viewed as a situation of ensembles inequivalence. 
Note, however, that it concerns metastable states (see Section 



interaction (see Section r2.4[) . Since the series of equilib- 
ria with positive magnetization does not present turning 
point of temperature, nor bifurcation, we conclude that 
the whole branch is stable. The series of equilibria with 
negative magnetization exists for t] > r]c{h). We know 
that the states with inverse temperature rj — > -|-oo and 
energy e'^mW ^^'^ stable since they arc local energy min- 
ima at fixed mass (see Section [2^ . On the other hand, the 
series of equilibria with negative magnetization presents 
a unique turning point of temperature at rjc{h). There- 
fore, starting from (ej„,;„(ft,), -t-co), the series of equilibria 
with negative magnetization is stable before the turning 
point of temperature CE and it becomes, and remains, 
unstable afterwards. By comparing the free energies of 
the stable solutions in competition (see Figure [5]), we see 
that the states with positive magnetization always have a 
lower free energy than the states with negative magneti- 
zation. Therefore, the states with positive magnetization 
(i.e. < a; < A-oo) are fully stable (S), the states with 
large negative magnetization (i.e. — oo < x < a;™"°(/i)) 
are metastable (M) and the states with small negative 
magnetization (i.e. a;™"°(/i) < x < XQ{h)) are unstable 
(U). In conclusion: (i) for rj < rjc{h) there exists a unique 
equilibrium state that is fully stable (S); (ii) for 77 > i]c{h) 
there exists one fully stable state (S), one metastable state 
(M) and one unstable state (U). 

Let us now consider the microcanonical ensemble. We 
have to determine the maxima of entropy at fixed mass 
and energy. A global maximum will be called fully sta- 
ble (S), a local maximum will be called metastable (M) 
and a minimum or saddle point will be called unstable 
(U). To apply the Poincare theorem, we have to plot 77 
as a function of e (sec Figure [T}. The series of equilibria 
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Fig. 11. Series of equilibria (caloric curve) giving the inverse 
temperature r; as a function of the energy e for different values 
of /i = 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. 




Fig. 12. Magnetization b (order parameter) as 
a function of the inverse temperature r] for h — 
0.001, 0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. 
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Fig. 13. Magnetization b (order parameter) as a function of 
the energy e for /i = 0.001,0.01,0.05,0.1,0.2,0.3,0.4,0.5,0.6. 
The dashed line corresponds to 6 = — (e/2)^''^. The "magic" 
points ±x,n where all the curves cross each other are defined 
in Section [5l 



with positive magnetization exists for any e > e,„i„(/i). 
We know that the system is stable at high energies since 
it becomes equivalent to a classical gas without interac- 
tion (see Section [531) • Since the series of equilibria with 
positive magnetization does not present turning point of 
energy, nor bifurcation, it follows that the whole branch 
is stable. The series of equilibria with negative magne- 
tization exists for emj„(/i) < e < (-c{h). We know that 
the states with energy e'^ini^) ^^'^ inverse temperature 
r] — >■ -t-oo are stable since they correspond to local en- 
ergy minima (see Section 12. 4p . On the other hand, the 
series of equilibria with negative magnetization presents a 
unique turning point of energy at edh). Therefore, start- 
ing from (e^j„(/i), +oo), the series of equilibria with neg- 
ative magnetization is stable before the turning point of 
energy MCE and it becomes, and remains, unstable af- 
terwards. By comparing the entropies of the stable solu- 
tions in competition (see Figure H]), we see that the states 
with positive magnetization always have a higher entropy 
than the states with negative magnetization. Therefore, 
the states with positive magnetization (i.e. < x < -l-oo) 
are full stable (S), the states with large negative mag- 
netization (i.e. — oo < a; < a;™*'^''°(/i)) are metastable 
(M) and the states with small negative magnetization (i.e. 
x'^^^^°{h) < X < xo{h)) are unstable (U). In conclusion: 
(i) for e > ec{h), there exists a unique equilibrium state 
that is fully stable (S); (ii) for eo(/i) < e < ec(/i), there 
exists one fully stable state (S), one metastable state (M) 
and one unstable state (U); (iii) for e'„j„(/i) < e < eo(/i), 
there exists one fully stable state (S) and one metastable 
state (M); (iv) for e„iin{h) < e < e^j„(/i), there exists 
only one fully stable state (S). 

If we only consider fully stable states (S), we conclude 
that the ensembles are equivalent and that the specific 
heat is always positive. Furthermore, there is no phase 
transition when h ^ 0, contrary to the case h = which 
displays a second order phase transition at (ec, Vc) = 
(1, 1). However, if we take into account metastable statefB 
we find a small region of ensembles inequivalence. Indeed, 
we note that the states situated in the region between 
points CE and MCE have negative specific heats (see Fig- 
ure [TU]). These states are unstable in the canonical ensem- 
ble while they are stable in the microcanonical ensemble. 
We recall that negative specific heat is a sufficient but not 
necessary condition of canonical instability. In particular, 
the states past point MCE are unstable (in both ensem- 
bles) while they have positive specific heats. We note that 
the series of equilibria becomes unstable in the canonical 
ensemble when the specific heat passes from positive to 
negative values (the point CE has infinite specific heat) 
while the series of equilibria becomes unstable in the mi- 
crocanonical ensemble when the specific heat passes from 
negative to positive values (the point MCE has zero spe- 
cific heat). 



Metastable states are very important in systems with long- 
range interactions because they have tremendously long life- 
times, scaling like . Therefore, they can be considered as 
stable states in practice [64165) . 
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Fig. 14. Zeroth order phase transitions corresponding to the 
reorganization of the system from the anti-aligned phase to the 
aligned phase. This is similar to the gravothermal catastrophe 
and isothermal collapse of self-gravitating systems [14] . 
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Fig. 15. Canonical phase diagram in the (h, rj) plane where h 
is the external parameter (magnetic field) and 77 is the control 
parameter (inverse temperature). 



These results are very similar to those obtained in the 
thermodynamics of self-gravitating systems (see [M] for 
a review). For self-gravitating systems, there also exists 
a region of negative specific heat in the microcanonical 
ensemble in which the temperature decreases when the 
energy increases. On the other hand, when the system 
reaches the turning point of temperature or energy (called 
spinodal points), the metastable branch disappears and 
the system undergoes an "isothermal collapse" (in the 
canonical ensemble) or a "gravothermal catastrophe" (in 
the microcanonical ensemble) until an equilibrium state is 
reached (for systems with a small-scale cut-off). Similar 
phenomena occur in the present problem (see Figure [Mj. 
In the canonical ensemble, when we reach the spinodal 
point ric{h), the metastable branch made of states with 
negative magnetization disappears and the system under- 
goes a sort of instability similar to the isothermal collapse 
(fixed T) . Similarly, in the microcanonical ensemble, when 
we reach the spinodal point tc{h), the metastable branch 
made of states with negative magnetization disappears 
and the system undergoes a sort of instability similar to 
the gravothermal catastrophe (fixed E). The system re- 
organizes itself and finally reaches an equilibrium state 
with positive magnetization. This corresponds to a zeroth 
order phase transition marked by the discontinuity of en- 
tropy and free energy (see Figs. |3][7]). There is no first or- 
der phase transition in the HMF model with a magnetic 
field contrary to the case of self-gravitating systems [2]. 
Furthermore, in the gravitational case, the system under- 
goes a transition from a homogeneous phase to a clus- 
tered phase while for the HMF model with a magnetic 
field, it undergoes a reorganization from an anti-aligned 
phase (magnetization pointing in a direction opposed to 
the field) to an aligned phase (magnetization pointing in 
the same direction as the field). 

Important remark: a more detailed stability analysis 
(see Appendix [D]) reveals that the metastable states arc, 
in fact, unstable with respect to perturbations that change 
the phase of the magnetization (i.e. 5By ^ 0). Therefore, 
if we allow for these perturbations, the anti-aligned phase 



becomes unstable and the richness of the problem disap- 
pears. In the following, we shall consider only perturba- 
tions for which 5 By = 0. However, one should keep in 
mind this remark in the interpretation of the results. 

2.7 Phase diagrams 

The preceding results can be summarized by drawing ap- 
propriate canonical and microcanonical phase diagrams 
(sec Figs, [m and [TH). 

Let us first describe the canonical phase diagram in 
the (/i, 77) plane represented in Figure [TS] The solid line 
corresponds to the critical inverse temperature ric{h) that 
exists only for h < 1 and tends to -l-oo when ft, — > 1^. 
This curve divides the parameter space in two regions. 
In the region denoted "aligned" , the magnetization has 
the same sign as the magnetic field. This is the case for 
any temperature when ft > 1 or for inverse temperatures 
rj < r]c{h) when ft < 1. The region denoted "mixed" cor- 
responds to a mixed zone in which the magnetization can 
be aligned or anti-aligned with the magnetic field. In the 
N ^ +00 limit, the system can be observed in only one 
of these two phases depending on the way it has been ini- 
tially prepared. For finite N, the system undergoes ran- 
dom transitions from the aligned phase (fully stable) to 
the anti-aligned phase (metastable) as described in Sec- 
tioning! As explained in Section all the stable states 
in the canonical ensemble have positive specific heat c > 0. 

The microcanonical phase diagram in the (ft, e) plane 
is represented in Figure 1161 The upper solid line corre- 
sponds to the critical energy ec(ft) that exists only for 
ft. < 1 and tends to 2 when ft — >■ 1~. We have also 
represented the characteristic energies eo(ft), ^'minW ^^^^ 
£minih). These curves divide the parameter space in two 
regions (of course, the region below the minimum energy 
^min[h) is forbidden). In the region denoted "aligned", the 
magnetization has the same sign as the magnetic field. 
This is the case for any energy e > eminih) when ft > 1 
and for energies e > ec{h) or emin{h) < e < ^'j^inW when 
ft < 1. The region denoted "mixed", delimited by the 
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Fig. 16. Microcanonical phase diagram in the (/i, e) plane 
where h is the external parameter (magnetic field) and e is 
the control parameter (energy). 



curves e'^i„{h) and edh), corresponds to a mixed zone 
in which the magnetization can be ahgned or anti-aligned 
with the magnetic field. In the N — !► +oo limit, the sys- 
tem can be observed in only one of these two phases 
depending on the way it has been initially prepared. 
The subpart of the mixed region delimited by the curves 
^minW ^o{h) exhibits phase space gaps responsible 



for ergodicity breaking as described in Section l4~2l There- 
fore, for finite N, the system undergoes random transi- 
tions between the aligned phase (fully stable) and the 
anti-aligned phase (metastable) when eo(/i) < e < edh) 
while it remains blocked in one of these two phases when 
^minW ^ £ < ^o{h) even for small N. On the other hand, 
the dashed line corresponds to the energy e{r]c{h)) associ- 
ated with the turning point of temperature. As explained 
in Section [2.5l the region between the curves e{r]c{h)) and 
ec{h) corresponds to states with negative specific heats 
c < which arc stable in the microcanonical ensemble 
but unstable in the canonical ensemble. This is therefore 
a region of ensembles inequivalence. 



3 Stability analysis in the canonical ensemble 

In this section, we analytically study the thcrmodynamical 
stability of the system in the canonical ensemble, using a 
procedure similar to the one developed in j61j for the case 
h = 0. 



3.1 The free energy F{B) 

The minimization problem (|T4)) determines the statisti- 
cal equilibrium state of the HMF model in the canonical 
ensemble. It is shown in Appendix A. 2. of [HB] that the 
solution of (|T4)) is given by 



1/2 



where p{6) is the solution of 

min {F[p] \ AI[p] 
p 



Af}, 



(42) 



(43) 



F[p] P^dG + J p^,,tde + T J plnpde, (44) 

is the configurational free energy. Therefore, the minimiza- 
tion problems ([T4| and (|43|) are equivalent: 



11 ^ (SI). 



(45) 



This equivalence holds for global and local minimization 
(i) f{6, v) is the global minimum of iff p{9) is the 
global minimum of (j43p and (ii) /(0, v) is a local minimum 
of (|14p iff p{9) is a local minimum of (|43p . We arc therefore 
led to considering the minimization problem (|43p which is 
simpler to study since it involves the density p{d) instead 
of the distribution function f(0, v). 

For the HMF model, the potential energy is given by 
equation (fTO)) so that the free energy (|44|) can be rewritten 



(46) 



ttB^ 27r 
F[p]^-—--B.H- 



T / p\npd0. 



Let us determine the global minimum of free energy at 
fixed mass. To that purpose, we shall reduce the minimiza- 
tion problem (|43p to an equivalent but simpler minimiza- 
tion problem. To solve the minimization problem (|43p , we 
proceed in two steps: wc first minimize F[p] at fixed M 
and Bx and By. Writing the variational principle as 

J p\np dO^ + aSM + n^SB^ + HySBy = 0, (47) 



wc obtain 



Pi{0)^Ae 



X COS 9 



(48) 



We have anticipated the fact that the magnetization of the 
global minimum of free energy is parallel to the magnetic 
field so that By = 0, implying py = 0. The Lagrange 
multipliers A — e^^^" and A = —^Px are determined by 
the constraints M and B^ (denoted B in the following) 
through the equations 



A = 



M 



27r/o(A) - 



and 



2ttB h{X) 



(49) 



(50) 



kM Io{X)' 

Equation (|48)) is the (unique) global minimum of F[p] with 
the previous constraints since S^F = J ^^j— dO > 
(the constraints are linear in p so that their second varia- 
tions vanish). Then, we can express the free energy F[p] as 
a function of B by writing F{B) = F[pi]. After straight- 
forward calculations, we obtain 



FiB) 



nB^ 



TX- 



2ttB 



MTlnloiX), 



(51) 
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where X{B) is given by equation (jSOp . FinaUy, the mini- 
mization problem (j43p is equivalent to the minimization 
problem 



mm{F(B)}. 



(52) 



in the sense that the solution of is given by equations 
(gnj) and dnO]) where B is the solution of Note 
that the mass constraint is taken into account implicitly 
in the variational problem ([5^ . Therefore, and (|52l) 
are equivalent for global minimization: 



13 ^ (E2 



(53) 



Furthermore, we show in Appendix ID. II that they are also 
equivalent for local minimization provided that we impose 
the constraint 6By = to the perturbations (otherwise 
the metastable states are always unstable). Under these 
conditions, the equivalence ([5^ holds for global and local 
minimization: (i) p(6) is the global minimum of (|43|) iff 
B is the global minimum of ([5^ and (ii) p{d) is a local 
minimum of ([^S]) iff i? is a local minimum of ([5^ . We 
are therefore led to considering the minimization problem 
(j52[) which is simpler to study since, for given T and M, we 
just have to determine the minimum of a function F{B) 
instead of the minimum of a functional F[p\ at fixed mass. 



3.2 The condition of canonical stability 

Let us therefore study the function F{B) defined by equa- 
tions ([FT]) and ([50]) for given T and M . Introducing the 
dimensionless variables of Section 12.21 we have to study 
the function 

/(6) = -26^ -46/1+ -A6- -ln/o(A), (54) 

where X{b) is given by equation (|50p and 77 is prescribed. 
Its first derivative is 

m = -4(6 + /.) + ^A + ^f5-^)§. (55) 
77 ?; V J ah 

Using the identity /o('^) = ^'iid the relation ([50]) . we 

see that the term in parenthesis vanishes. Then, we get 



/'(&) = -4(6 + /i) + -A. 



(56) 



The critical points of /(6), satisfying /'(6) = 0, correspond 
therefore to 



\^ X = 2ri{b + h). 



(57) 



Substituting this result in equation (|50p . we obtain the 
self-consistency relation 



h{2i^{b + h)) 

h{2i^{b + h)y 



(58) 



which determines the magnetization 6 as a function of the 
inverse temperature rj. This returns the equilibrium results 
of Section O 
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Fig. 17. Graphical construction determining the critical points 
of /(t() and their stability. The critical points are determined 
by the intersection(s) between the curve h — 6(A) defined by 
equation (|50p and the straight line b — X/{2i]) — h. The critical 
point is a minimum (resp. maximum) of f{b) if the slope of the 
curve 6(A) at that point is smaller (resp. larger) than the slope 
of the straight line 6 = \/{2r^) — h. 



Now, a critical point of f{b) is a minimum if f"{b) > 
and a maximum if f"(b) < 0. Differentiating equation (|56p 
with respect to 6, we find that 

/"W = -§-4. (59) 
r] ab 

Therefore, a critical point A = a; is a minimum if 

b\x) < i-, (60) 

and a maximum if the inequality is reversed. According to 
equation (j37p . the turning point of temperature (drj/dx = 
0) corresponds to 



b'ix) 



1 



2r]{x) 



(61) 



Comparing equation (|5T|) with equation (pO)) . we see that 
the change of stability in the canonical ensemble corre- 
sponds to the turning point of temperature ric{h) occuring 
at x^^^°{h) in agreement with the Poincare theorem. Us- 
ing equation (|4T|) . it is also easy to establish that when the 
specific heat is negative, the converse of inequality is 
always fulfilled so that the system is unstable in agreement 
with general theorems of statistical mechanics. 

We can determine the minima and maxima of the func- 
tion f{b) by a simple graphical construction. To that pur- 
pose, we plot b 0jS 0, function of A according to 6(A) = 
/i(A)//o(A). This function is represented in Figure [T71 We 
note that 6 — !> ±1 for A — > ±00 and that 6 ^ A/2 for 
A ^ 0. According to equation (j57[) . the critical points of 
/(6) are determined by the intersection of this curve with 
the straight line 6 — \/{2r]) ~ h. For given 77, this deter- 
mines x{ri) and 6(77). For 77 < r]c{h), there is a unique so- 
lution 65 > that has the same sign as the imposed mag- 
netic field. For j] — 0, we see that xs 0+ and 65 — > 0+. 
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Fig. 18. Free energy f{b) as a function of the magnetization 
b for a given value of the inverse temperature rj and magnetic 
field h (specifically h = 0.2). For rj < r)c(h) ~ 1.89, this curve 
has a (unique) global minimum at bs > 0. For 77 > r]c{h), this 
curve has a global minimum at fes > 0, a local minimum at 
6m < and a local maximum at &[/ < 0. 



For 77 > ric{h), there are three solutions: a solution 65 > 
that has the same sign as the imposed magnetic field and 
two solutions bu < and bjxi < bu whose sign is opposite 
to the sign of the imposed magnetic field. For 77 — > +00, 
we see that xs +00, 65 — ?> 1, x^/ — > —00, b^j — > — 1, 
xu xo{h) and bjj — > bo{h) = —h. According to inequal- 
ity (pO)) . a solution is a minimum of /(6) if b'{x) < 1/(277) 
and a maximum if b'{x) > 1/(277). Therefore, a critical 
point of free energy f(b) is a minimum (resp. maximum) 
if the slope of the main curve is lower (resp. higher) than 
the slope of the straight line at the point of intersection. 
From this criterion, we directly conclude that the solutions 
bs and 5a/ are minima of free energy while the solution 
bu is a maximum of free energy. 

To complete our analysis, it can be useful to plot the 
function f{b) for a prescribed inverse temperature rj. It is 
given in parametric form by 

/(A) = -2b{\f - 4:b(\)h+-Xb(X) - -ln/o(A), (62) 

?7 77 

together with equation (|50p . Eliminating A between the 
expressions ([5^ and ([50)) . we obtain the free energy f{b) 
as a function of the magnetization b for a fixed value of 
the inverse temperature 77. For rj < r/c{h) and 77 > rjc{h), 
this function displays the two behaviors described above, 
as illustrated in Figure [TSl 

It can be shown that the probability of a fluctuation 
with magnetization 6 in the canonical ensemble is given by 
PcE{b) = z(;^e-^''^-^('') (see Appendix E]) . When there 
exists only one (global) minimum of free energy bs^ as in 
the case /i > 1 or in the case h <1 and 77 < ric{h), the situ- 
ation is simple. For N — >■ +00, the distribution of magneti- 
zation is strongly peaked around the optimal value b = bs 
so that, after a transient regime, the system will be found 
in that state. For h <1 and 77 > r]c{h), the free energy f[b) 
presents a global minimum at 65 and a local minimum at 
For -t-oo, these states have infinite lifetime and 



the system will be found in one of them depending on how 
it has been initially prepared. For finite N , the system will 
jump from one state to the other. Of course, it will spend 
more time in the global minimum of free energy (fully sta- 
ble) than in the local one (metastable). To pass from one 
state to the other, the system has to overcome an entropy 
barrier which is played by the unstable solution bjj. This 
barrier scales like ^iiN\Af\. For finite A, this barrier is fi- 
nite so that random transitions from the fully stable state 
bs to the metastable state 6a/ are possible. For A^ -|-oo, 
the barrier is too hard to cross and the system remains in 
one of these two states as previously indicated (according 
to the Kramers formula, the lifetime of these states scales 
like tufe - e^"^l^^l ^ +00) [MI55] . 

4 Stability analysis in the microcanonical 
ensemble 

In this section, we analytically study the thermodynamical 
stability of the system in the microcanonical ensemble, 
using a procedure similar to the one developed in [61) for 
the case /i = 0. 



4.1 The entropy S{B) 

The maximization problem ([T3| determines the statistical 
equilibrium state of the HMF model in the microcanonical 
ensemble. It is shown in Appendix A.l. of that the 
solution of ((T3)) is given by 

fie,v) = (^-l-y\{e)e-^'^, (63) 

where the inverse temperature /? = 1/T is determined by 
the energy constraint 

E = ^MT + W, (64) 

and p{6) is the solution of 

msix{S[p]\M[p]^ M}, (65) 
p 

where 

S[p] ^^MlnT- J p\npd0, (66) 

is the configurational entropy. Eliminating the tempera- 
ture thanks to the constraint ()64[). we can write the en- 
tropy in terms of p alone as 

S[p] = -J plnpde + ^M\n{E ~W[p]). (67) 

Therefore, the maximization problems p3p and (|65l) are 
equivalent: 

(dSl) <;=^ (ESI). (68) 
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This equivalence holds for global and local maximization 
[66] : (i) f{d, v) is the global maximmii of (fT3|) iff p{9) is the 
global maximum of (j65p and (ii) f{6,v) is a local maxi- 
mum of (fT51) iff p{6) is a local maximum of We arc led 
therefore to considering the maximization problem (|65p 
which is simpler to study since it involves the density p{9) 
instead of the distribution function v). 

For the HMF model, the potential energy is given by 
equation ([TO)) so that the energy ([M|) and the entropy ([57]) 
can be rewritten 



2 



(69) 



S[p] = -J plnpde+^AIln (^E 



2n 



(70) 



Let us determine the global maximum of entropy at fixed 
mass (the conservation of energy is automatically taken 
into account in equation (j70|) '). To that purpose, we shall 
reduce the maximization problem (j65p to an equivalent 
but simpler maximization problem. To solve the maxi- 
mization problem (j65p . we proceed in two steps: we first 
maximize S[p] at fixed M and B^ and By. Writing the 
variational problem as 

-^(/„„..,)-„™ 0,(71, 

we obtain 



Pi{e)^Ae 



A cos 6 



(72) 



where A and A are determined by the constraints M and 
B through the equations 



M 



and 



b = 



2TrB _ /i(A) 
fcM ~ /o(A)' 



(73) 



(74) 



We have anticipated the fact that the magnetization of 
the global maximum of entropy is parallel to the magnetic 
field so that = 0, implying py = 0. Equation ([7^ is 
the (unique) global maximum of S[p] with the previous 

constraints since S^S = — 5 / ^^''p d6 < (the constraints 
are linear in p so that their second order variations vanish). 
Then, we can express the entropy 5 as a function of B by 
writing S{B) = S[pi]. After straightforward calculations, 
we obtain 



SiB)=MlnIo{X) - ^A 



2 



E 



ttB^ 



k 



(75) 



where X{B) is given by equation (|74p . Finally, the maxi- 
mization problem (j65p is equivalent to the maximization 
problem 



max{5(i?)} 



(76) 



in the sense that the solution of (|65|) is given by equa- 
tions ([72])- (HI) where B is the solution of ((75]). Note that 
the energy and mass constraints arc taken into account 
implicitly in the variational problem (j76p . Therefore, (|65p 
and (j76p are equivalent for global maximization: 



§5\) ^ (HH). 



(77) 



Furthermore, we show in Appendix ID . 2 1 that they are also 
equivalent for local maximization provided that we impose 
the constraint SBy = to the perturbations (otherwise 
the metastable states are always unstable). Under these 
conditions, the equivalence ([77]) holds for global and local 
maximization: (i) p{6) is the global maximum of (|65p iff 
B is the global maximum of ([76)) and (ii) p{9) is a local 
maximum of ([55)) iff S is a local maximum of ([75)) . We 
arc therefore led to considering the maximization problem 
([76)) which is simpler to study since, for given E and M, we 
just have to determine the maximum of a function S{B) 
instead of the maximum of a functional S[p] at fixed mass. 



4.2 The condition of microcanonical stability 

Let us therefore study the function S{B) defined by equa- 
tions ()75p and ()74p . Introducing the dimensionless vari- 
ables of Section [2?2l we have to study the function 

s{b) = ln/o(A) - 6A + i In (e + 2b^ + ibh) , (78) 

where A(6) is given by equation ()74l) . Its first derivative is 

2{b+h) 



,,,, //^(A) \ dX , 
'^'^=[-m-Vd-b-'+e + 2b-^ + Abh 



, (79) 



Using the identity Io{x) = Ii{x) and the relation ([71)) . we 
see that the term in parenthesis vanishes. Then, we get 



s'ib) = ~X + 2ri{b + h), 



(80) 



where the inverse temperature 77 is determined by the en- 
ergy constraint ()69p which can be rewritten in dimension- 
less form 



€=--2b^~ Ahb. 



(81) 



The critical points of s(5), satisfying s'{b) = 0, correspond 
to 



A = a; = 277(6 + /i). 



(82) 



Substituting this result in equation (I74L we obtain the 
self-consistency relation 



Jl(27?(&+/l)) 

Ioi2rib + h))' 



(83) 
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which, together with equation (j8ip . determines the mag- 
netization as a function of the energy. This returns the 
equihbrium relationships of Section 12.21 

Now, a critical point of s{b) is a maximum if s"{b) < 
and a minimum if s"{b) > 0. Differentiating equation (|80p 
with respect to b, and recalling that the inverse temper- 
ature 77 is a function of b given by equation ()81|) . we find 
that 

s"{b) = 2^-^-8v'{b + hf. (84) 
Therefore, a critical point of s(6) is a maximum if 

^ > 277 -877^(6 + (85) 



b'{x) 



and a minimum if the inequality is reversed. Using equa- 
tion (15^ . inequality ([55)1 can be rewritten 



1 



b'{x) 



> 2(77 - a;2). 



(86) 



According to equation (|38p . the turning point of energy 
{de/dx = 0) corresponds to 



6' (a;) 



(87) 



Comparing equation ((87|) with equation (|86)) . we see that 
the change of stability in the microcanonical ensemble cor- 
responds to the turning point of energy ec{h) in agreement 
with the Poincare theorem. On the other hand, using equa- 
tion ([5^. we note that 
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Fig. 19. Accessible ranges of magnetizations when h < 1 
(specifically h = 0.2). 
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Fig. 20. Accessible ranges of magnetizations when h > 1 
(specifically h = 2). 



s"{b) 



8if(b + hf 



Since the last term in equation (|88[) is negative, we recover 
the fact that canonical stability implies microcanonical 
stability . Indeed, if the critical point is a minimum of 
free energy {f"{b) > 0), then it is a fortiori a maximum 
of entropy (s"(6) < 0). 

To complete our analysis, it can be useful to plot the 
function s(6) for a prescribed energy e. It is given in para- 
metric form by 

s(A) = ln/o(A) - 6(A)A + ^ In (e + 2b{X)^ + Ahb{X)) , 

(89) 

together with equation (|74l) . Eliminating A between equa- 
tions and (j74p . we obtain the entropy s{b) as a func- 
tion of the magnetization b for a fixed value of the energy 
e. Before going further, we must take into account the pos- 
sibility of phase space gaps in the system. Indeed, since 
the temperature is positive, the energy equation (j8ip im- 
plies that e + 2b^ + 4:hb > 0. Depending on the value of the 
energy, this constraint may restrict the range of accessible 
magnetizations (see Figs. [TOl and EO)) . We must distinguish 
several cases. 



Let us first assume that h < 1. (i) For e > edh), 
all magnetizations are allowed and the entropy s{b) has 
a unique (global) maximum at 65 > 0. (ii) For eo{h) < 
e < ec(h), all magnetizations are allowed and the entropy 
s{b) has a global maximum at 6s > 0, a local maximum 
at 6a/ < and a local minimum &i bu < 0. (iii) For 
e'minW < £ < co(^)i only magnetizations in the intervals 
[-l,6_(/i)[ and ]6+(/i), 1] with 



b±{h,e) 



-h±hjl 



(90) 



are allowed and the entropy s(6) has a global maximum 
at 6s > and a local maximum at 6j\/ < 0. (iv) For 
^minih) < e < e'minWj o^iIy magnetizations in the range 
]6+(/i),l] are allowed and the entropy s(6) has a unique 
(global) maximum at 6s > 0. These results are illustrated 
in Figure [5T] 

Let us now assume that h > 1. {i) For e > i'minW^ 
all magnetizations are allowed and the entropy s(6) has a 
unique (global) maximum at 6s > 0. (ii) For trnin{h) < 
e < e'riiinih), only magnetizations in the range ]6+(/i),l] 
are allowed and the entropy s(6) has a global maximum 
at 6s > 0. These results are illustrated in Figure 
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Fig. 21. Entropy s{b) as a function of the magnetization for 
a given value of the energy e when h < 1 (specifically h = 
0.2). For e > ec{h) ~ 0.458, this curve has a unique (global) 
maximum at bs > 0. For eo{h) = 0.08 < e < ^c{h), it has a 
global maximum at &s > 0, a local maximum at 6a/ < and a 
local minimum at bu < 0. For trninW = —1.2 < e < to{h), it 
has a global ma:ximum at 6s > and a local maximum at bM < 
0. For tmin{h) — —2.8 < e < e'mi„{h), it has a unique (global) 
maximum at bs > 0. For e'„^in{h) < e < eo{h), there is a gap of 
magnetization corresponding to the interval [6_ {h, e), b+{h, e)] 
and for emin{h) < t < e^i„(/i), only the interval ]b+{h, e), 1] is 
accessible. 
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b 

Fig. 22. Entropy s(6) as a function of the magnetization for 
a given value of the energy e when h > 1 (specifically h = 
2). This curve has a unique (global) maximum at bs > 0. 
For tmin{h) = — 10 < e < t'minih) = 6, only the interval 
]b+{h, e), 1] is accessible. 



+00, these states have infinite lifetime and the sys- 
tem will be found in one of them depending on how it has 
been initially prepared. For finite N, the system wiU jmnp 
from one state to the other. Of course, it will spend more 
time in the global entropy maximum (fully stable) than 
in the local one (metastable). To pass from one state to 
the other, the system has to overcome an entropy barrier 
which is played by the unstable solution bir. This bar- 
rier scales like 7V|zis|. For finite N, this barrier is finite 
so that random transitions from the fully stable state bs 
to the metastable state 6a/ are possible. For N — >■ -|-oo, 
the barrier is too hard to cross and the system remains 
in one of these two states as previously indicated (accord- 
ing to the Kramers formula, the lifetime of these states 
scales like t^/e ^ e^^* — -l-oo) |64l65j . Interestingly, for 
^'minW < ^ < ^o{h), the entropy curve still possesses 
two entropy maxima bs and bj^i but the model exhibits 
a gap in the magnetization corresponding to the inter- 
val [b-{h, e),b+{h, e)]. Such phase space gaps have been 
previously observed in other systems with long range in- 
teraction and they are responsible for ergodicity breaking 
[68l69j . Indeed, in the presence of such a gap, the system 
cannot jump from one state to the other even when N 
is small (the gap is equivalent to having an infinite en- 
tropy barrier between the two states). In that case, the 
system remains blocked in one of these two states for 
ever and ergodicity is broken (time averages do not co- 
incide with ensemble averages). This type of ergodicity 
breaking has been illustrated in [68169] for a generalized 
isotropic XY model with two and four mean-field inter- 
actions. The HMF model with a magnetic field is another 
system (maybe simpler) where ergodicity breaking should 
be observed in some range of parameters. Note, in con- 
trast, that these features do not arise in the canonical 
ensemble (see Section 13. 2p since all the values of the mag- 
netization are accessible. 

Remark: in the usual HMF model {h = 0), the entropy 
maximum is degenerate due to the rotational U{1) sym- 
metry. Therefore, there exists an infinity of equilibrium 
states that only differ by their phase (j) or equivalently by 
the position of the density maximum. For finite N, the 
system will explore these different maxima randomly. For 
N — > +00, it will remain blocked in one of them. 



5 Magnetic susceptibility 



It can be shown that the probability of a fiuctuation 
with magnetization b in the microcanonical ensemble is 
given by PMCE{h) = (sec Appendix |F|. When 

there exists only one (global) entropy maximum bs, as 
in the case ft, > 1 or in the case h < 1 and e > ec{h) 
or emin{h) < e < KninW- situation is simple. For 
N ^ +0O, the distribution of magnetization is strongly 
peaked around the optimal value 6 = 6s so that, af- 
ter a transient regime, the system will be found in that 
state. For h < 1 and eo{h) < e < £c{h), the entropy 
curve presents a global maximum at 6s and a local max- 
imum at 6j\/ separated by a local minimum ar bjj- For 



5.1 Canonical ensemble 

The curve giving the magnetization 6 as a function of the 
magnetic field h at fixed inverse temperature rj (canonical 
ensemble) is given in parametric form by 



2ttB _ Ii{x) 
Jm ~ Io{x)' 



2tiH _ X 
kM ~ 2r] 



b{x) 



(91) 



(92) 
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The magnetic susceptibiHty in the canonical ensemble is 
defined by x = {dB / dH)T- Introducing dimensionless 
variables, it can be written 



dhj, 



Using equations (|9T|) and ((92)) . we obtain 



2r, 



(93) 



(94) 



where b'{x) is given by equation We observe that 

the magnetic susceptibility (|94)) diverges when the condi- 
tion (j6ip is fulfilled, i.e. at the turning point of tempera- 
ture ric{h). Furthermore, we have the following results: (i) 
the stable branch (S) has positive magnetic susceptibility 
for all a; > 0, i.e. for all inverse temperatures rj. (ii) The 
metastable branch (M) has positive magnetic susceptibil- 
ity for all X < a:™"°(/i), i.e. for all 77 > r]c{h). (iii) The 
unstable branch (U) has negative magnetic susceptibility 
for all x™"°(/i) < X < xo{h), i.e. for all rj > ?]c{h). These 
results are illustrated on Figure [531 

It is a general result of statistical mechanics that the 
magnetic susceptibility is positive in the canonical ensem- 
ble (x > 0) since it is a measure of the variance of the 
fluctuations of the magnetization x = ^/5((^^)^)- Simi- 
larly, the specific heat is positive in the canonical ensemble 
since it is a measure of the variance of the fiuctuations of 
energy C = (3'^{{AEY). For the HMF model with a mag- 
netic field, our study shows that, in the canonical ensem- 
ble, all unstable states (a;^""°(/i) < x < xo{h)) have neg- 
ative magnetic susceptibility but only a fraction of them 
{x'^°-^°{h) < X < x^^''^°{h)) has negative specific heats 
(see Section 153)) . 

We shall now study the magnetic curve b{h) for a fixed 
inverse temperature rj. To understand the following dis- 
cussion, it can be useful to consider Figure [1^ in parallel. 
On this figure, we fix the value of 77 and progressively in- 
crease the value of h, starting from h = 0. Let us introduce 
the critical magnetic field hcij]) such that r]cihc) = V- Foi' 
■q < Tic = 1, this equation has no solution. For rj > rjc = I, 
this equation has one solution. We are now ready to study 
the curves b{h) for different values of the inverse temper- 
ature 77. 

li rj < rjc = 1: there is only one stable state 65 > 
with xs > 0. 

If 77 > ?/c = 1: (i) For < h < /ic(?/); there is one stable 
state bs > with xs > 0, one metastable state 6m < 
with xm > and one unstable state bu < with xu < 0- 
(ii) For h > hc{'q), there is only one stable state bs > 
with Xs > 0. 

Of course, the magnetic curves b{h) are antisymmetric 
with respect to a change of sign h —h oi the magnetic 
field. They are represented in Figures [M] and [25] for the 
different cases described above. For rj < rjc = 1 (see Fig- 
ure I24p . the curve is uni valued, going from 6 = — 1 for 
h -> —00 to 6 = -1-1 for h — >• +00. For rj > rjc = 1 (see Fig- 
ure [221), the curve b{h) displays an hysteretic cycle similar 
to the one observed for the Ising model in a magnetic field. 
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Fig. 23. Magnetic susceptibility x versus inverse temperature 
r\ for a magnetic field h = 0.01. Using the results of Appendix 
1X1 we find that x '7 for — >■ (corresponding to a: — >■ 0), 
X ~ 1/(4(1 -|- /i)^??) for 77 — >■ +CXJ (corresponding to a:: — !> -Hoo), 
X ~ 1/(4(1 — ^)^'7) for 77 — i- -l-oo (corresponding to a; — >■ —00) 
and X ~1 for t\ — >■ -l-oo (corresponding to x ^ xo{h)). On 
the other hand, x oc iiz{rj — rjc{h))~^^^ when 77 — >■ r]c{h)^ (corre- 
sponding to a; ^- x''^""°{h)). The stable (S) and metastable (M) 
states have x > while the unstable (U) states have x < 0. 
Note that the susceptibility increases rapidly close to ric{h). 
This is related to the second order phase transition at r/c = 1 
when h = Q (see Section [5. 3p . 
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Fig. 24. Magnetization h versus magnetic field h in the canon- 
ical ensemble for 77 < t^c = 1. 



For large positive h, the magnetization tends to its max- 
imum value 6 = -|-1. When we reduce the magnetic field 
h. the magnetization decreases since x > 0- For h > 0, 
the upper branch (S) is fully stable and corresponds to 
states whose magnetization 6 > has the same direction 
as the magnetic field. When /i < 0, this branch becomes 
metastable (M) since these states have a magnetization 
opposite to the magnetic field. However, since metastable 
states have tremendously long lifetimes for long-range in- 
teractions j64l65j . the system is expected to remain on 
this branch. For h < —hc{r]), the metastable branch dis- 
appears and the system jumps on the fully stable branch 
(S) with a magnetization 6 < having the same direction 
as the magnetic field. This transition corresponds to the 
reorganization of the system from an anti-aligned phase 
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Fig. 25. Magnetization h versus magnetic field h in the canon- 
ical ensemble for r\ > r\c = \. This curve displays a classical 
hysteretic cycle similar to the Ising model in a magnetic field. 



to an aligned phase. If we keep decreasing the magnetic 
field, the magnetization decreases since x > until the 
minimum value b ~ —1 reached for h — >■ — oo. If we now 
increase the magnetic field, the magnetization increases. 
For h < 0, the system is on the stable branch (S) cor- 
responding to states with negative magnetization 6 < 
aligned with the magnetic field. For h > 0, this branch 
becomes metastable (M) since these states have a mag- 
netization opposite to the magnetic field. However, the 
system remains on this branch for the reason given pre- 
viously. For h > hcij]) the metastable branch disappears 
and the system jumps on the fully stable branch with pos- 
itive magnetization. We have thus described a classical 
hysteretic cycle. In the canonical ensemble, the stable and 
metastable states have positive magnetic susceptibilities 
and unstable states have negative susceptibilities like for 
the classical Ising model. 



5.2 Microcanonical ensemble 

The curve giving the magnetization 6 as a function of the 
magnetic field h at fixed energy e (microcanonical ensem- 
ble) is given in parametric form by 



h{x) 



_ 27rB _ h{x) 
e - f 6(x) + 26(x)^ 



(95) 



(96) 



The magnetic susceptibility in the microcanonical ensem- 
ble is defined by % = [dB/dH)E- Introducing dimension- 
less variables, it can be written 



X 



[dh 



Using equations (pSj) . and ([5^ . we obtain 
2t]{x) (1 - 2b{x)x) 



Xix) 



(97) 



(98) 



Fig. 26. Magnetic susceptibility x versus energy e for a mag- 
netic field h = 0.01 < hm- We find that x ~ 1/e for £ ^ +oo 
(corresponding to a; ^ 0), x ^ ~1/(2(1 + h)) for e — > eminih) 
(corresponding to a; — i- +oo), x ^1/(2(1 — h)) for e — t- 
^'min{h) (corresponding to a- — s> —oo) and x ~^ ~1 — 2hxo{h) 
for e eo(ft) (corresponding to x — i- xo{h)). On the other 
hand, x oc ±(ec(ft) — e)~^^^ when e — >■ ec(h)~ (correspond- 
ing to a; — > x^^"'° (h)) . The stable and metastable states have 
X > for e > Em and x < for e < e^. Note that the suscep- 
tibility increases rapidly close to ec{h). This is related to the 
existence of the second order phase transition at = 1 when 
h — (see Section r5.3|l . 



where b'{x) is given by equation pop. In the microcanoni- 
cal ensemble, the magnetic susceptibility x of stable states 
can be either positive or negative [1]. Indeed, it can be 
shown that x = ^{I3{{AB)^) + B{dB/dE)) [52 (see also 
Appendix IF|) . We first observe that the magnetic suscep- 
tibility diverges when the condition ([57| is fulfilled, i.e. at 
the turning point of energy. We also observe that the mag- 
netic susceptibility vanishes when b{x) = 1/ (2x). Together 
with equation ([95]) this gives x = ±Xm — ±1.0657. Then, 
using equations ([5T|l - ([55|) . we obtain b = ±6^ — ±0.46918 
and e = Cm = l/(2a;^„) — 0.44025. We stress that these 
values are independent ofh! (by contrast, the correspond- 
ing inverse temperature rj = ?7,„(/i) ~ 0.53285/(0.46918 + 
h) depends on h). Therefore, all the series of equilibria 
6(e) with positive magnetization pass by the "magic" point 
(em, bm)- On the other hand, all the series of equilibria 6(e) 
with negative magnetization pass by the "magic" point 
(cm,— 6m) provided that < h < bm (corresponding to 
Xoih) > —Xm)- Let us finally determine the relative posi- 
tion of the turning point of energy x™'^'^'^°{h) with respect 
to the point —Xm at which the magnetic susceptibility 
vanishes. A simple calculation shows that they coincide 
for h = h,n = l/{x„ii4:Xm - 1)) - 0.2649. More precisely, 
x™*'=™(/i) > -x,n for 0<h<hm while x^'^'^^ih) < -Xm 
for h,n < h < bm- These results can be checked on Figure 

m 

We can now state the main results (it is recommended 
to consider Figure [Ql in parallel): (i) for any h > 0, the 
stable branch (S) has positive magnetic susceptibility for 
< a; < Xm (i.e. e > e™) and negative magnetic sus- 
ceptibility for X > Xm (i.e. eminih) < e < e„,); (ii) for 
h < hm, the metastable branch (M) has negative magnetic 
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Fig. 27. Magnetic susceptibility x versus energy e for a mag- 
netic field hm < h = O.Z < hm- The metastable states have 
X < 0. The stable states have x > for e > and x < for 
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Fig. 29. Magnetization h versus magnetic field h in the micro- 
canonical ensemble for ec = 1 < e = 1.25 < 2. 
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Fig. 28. Magnetic susceptibility x versus energy e for a mag- 
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The stable states have x > for e > and x < for e < 6^- 



susceptibility for x < —Xm (i.e. ^'minW ^ ^ ^ ^m) and 
positive magnetic susceptibility for —Xm < x < 
(i.e. Em < e < ^c{h)). The unstable branch (U) has neg- 
ative magnetic susceptibility for 2;™*'^''°(/i) < a; < XQ{h) 
(i.e. eo{h) < e < ^c{h)); (iii) for < h < hm, the 
metastable branch (M) has negative magnetic susceptibil- 
ity for X < x^'^'-^ih) (i.e. e'^^^{h) < e < edh)). The un- 
stable branch (U) has positive magnetic susceptibility for 
^micro^^^ < a; < —Xm (i-G. £,« < e < (c{h)) and negative 
magnetic susceptibility for —Xm < x < xo{h) (i.e. eo{h) < 
e < Cm); (iv) toT bjn < h < 1, the metastable branch (M) 
has negative magnetic susceptibility for x < x™'^'^°{h) (i.e. 
e'minih) < e < f-c{h))- The unstable branch (U) has pos- 
itive magnetic susceptibility for a:J?"^''°(/i) < a; < XQ{h) 
(i.e. eo(/i) < e < f-c[h)). These results are illustrated on 
Figures HHi 

We shall now study the magnetic curve h{h) for a fixed 
energy e. To understand the following discussion, it can be 
useful to consider Figure [T51 in parallel. On this figure, we 
fix the value of e and progressively increases the value of 
h, starting from /i = 0. To prepare the following discus- 
sion, we introduce some critical magnetic fields. Let /ic(e) 



be the magnetic field such that €c{hc) = e. For e > 2, this 
equation has no solution. For ec = 1 < e < 2, this equation 
has a unique solution. For e,„ < e < ec = 1, this equation 
has two solutions denoted hi^\e) and h'^^e). For e < £„, 
this equation has no solution. Let /io(e) be the magnetic 
field such that eg {ho) = e. For e > 2, this equation has 
no solution. For < e < 2, this equation has a unique 
solution explicitly given by ft.o(e) = (e/2)^^^- For e < 
it has no solution. Let /im,j„(e) be the magnetic field such 
that i'mini^'min) = ^- ^'^^ £ > 2, this cquation has no solu- 
tion. For — 2 < e < 2, this equation has a unique solution 
explicitly given by /i^„;„(e) = (e + 2)/4. For e < —2, this 
equation has no solution. Let /irnMi(e) be the magnetic field 
such that e-min{hmin) = £• For e < —2, this equation has a 
unique solution explicitly given by /imm(e) = —{(- + 2)/4. 
For e > —2, this equation has no solution. We are now 
ready to study the curves b{h) for different values of en- 
ergy e. 

If e > 2: there is only one stable state 6s > with 
Xs > 0. The magnetic curve b{h) is similar to the one of 
Figure [M] in the canonical ensemble. 

If ec = 1< e < 2 (see Figure[Ml): (i) For < /i< hc{e), 
there is only one stable state 65 > with xs > 0. (ii) For 
/ic(e) < h < ho{e), there is one stable state bs > with 
Xs > 0, one metastable state 6a/ < with xm < and one 
unstable state bjj < with xu > 0. (iii) For ho{e) < h < 
/i^j„(e), there is one stable state 65 > with xs > and 
one metastable state 6m < with xm < 0. (iv) For h > 
h'mini^)' there is only one stable state 65 > with xs > 0. 
If we start from the stable branch (S), we remain on this 
branch for all h. Furthermore, the magnetization increases 
with h since x > 0. If we start from the metastable branch 
(M), made of anti-aligned states, then we jump on the 
stable branch, made of aligned states, for \h\ > h'^^^{e) or 
\h\ < hc{e). Note that the metastable branch has negative 
susceptibility so that 6 decreases with h. 

If e,n < e < Ec = 1 (see Figure 15(11) : (i) For < h < 
hi^\€), there is one stable state 65 > with xs > 0, one 
metastable state 6j\/ < with xm > and one unstable 
state bu <0 with xu < 0. (h) For h[^\e) < h < hl?\e), 
there is only one stable state 65 > with xs > 0. (iii) 
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Fig. 30. Magnetization b versus magnetic field h in the micro- Fig. 32. Magnetization b versus magnetic field h in the micro- 
canonical ensemble for em < t ~ 0.5 < ec — 1. canonical ensemble for —2 < e = —1 < 0. 
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For h)r'{e) < h < ho{e), there is one stable state 65 > 
with xs > 0, one metastable state bM < with xm < 
and one unstable state bu < with xu > 0. (iv) For 
ho{e) < h < h'^^^{e), there is one stable state 65 > 
with Xs > and one metastable state 6m < with xm < 
0. (iv) For h > h'^^^{e), there is only one stable state 
bs > with xs > 0. The central curve displays a classical 
hysteretic cycle similar to the one described in Section 
15.11 in the canonical ensemble. On the other hand, if we 
start from the metastable branch on the periphery, then 
the system undergoes a transition to the stable branch for 

\h\ > h'^^^{e) or \h\ < he (e) as in the previous case. 

If < e < e„ (see Figure EH): (i) For < /i < ho{e), 
there is one stable state bs > with xs < 0, one 
metastable state 6a/ < with xm < and one unstable 
state 6(7 < with xu < 0. (ii) For ft,o(e) < h < /i^i„(e), 
there is one stable state 65 > with xs < and 
one metastable state 6m < with xm < 0. (iii) For 
h > /imi„(e), there is only one stable state 6s > with 
Xs < O.This case is interesting since it yields a new type 
of hysteretic cycle related to the fact that the stable and 
metastable states have negative magnetic susceptibility. 
For large positive h, the magnetization tends to 6 = 6s. 
When we reduce the magnetic field h, the magnetization 



increases since x < 0- For ^ > 0, the upper branch (S) 
is fully stable and its magnetization 6 > has the same 
direction as the magnetic field. When h < 0, this branch 
becomes metastable (M) since its magnetization is now 
opposite to the magnetic field. However, since metastable 
states have tremendously long lifetimes for long-range in- 
teractions f64l65j , the system remains on this branch. The 
magnetization increases until the maximum value 6 = 1 so 
that the system is in its most anti-aligned configuration. 
For h < — /i!„i„(e), the metastable branch disappears and 
the system jumps on the fully stable branch (S) with a 
magnetization 6 < having the same direction as the 
magnetic field. This transition corresponds to the reorga- 
nization of the system from an anti-aligned phase to an 
aligned phase. If we decrease the magnetic field, the mag- 
netization increases to — 6^ since x < 0. If we now increase 
the magnetic field, the magnetization decreases. For h < 0, 
the system is on the stable branch (S) corresponding to 
states with negative magnetization 6 < aligned with the 
magnetic field. For h > 0, this branch becomes metastable 
(M) since these states become anti-aligned with the mag- 
netic field. However, the system remains on this branch 
for the reason given previously. For h > h'^^^{e) the 
metastable branch disappears and the system jumps on 
the fully stable branch with positive magnetization. We 
have thus described a new type of hysteretic cycle. In the 
microcanonical ensemble, the stable and metastable states 
have negative magnetic susceptibilities so that, increasing 
the magnetic field, the magnetization decreases. 

If -2 < e < (see Figure[32l): (i) For < /i< h'^^Je), 
there is one stable state 6s > with xs < 3'iid 
one metastable state 6m < with xm < 0. (ii) For 
h > /imm(c)j there is only one stable state 6s > with 
Xs < 0. The discussion is essentially the same as the one 
given previously since the only difference is the nonexis- 
tence of unstable states (which play no role) in the present 
case. 

If e < -2 (see Figure (i) For < /i < /i™„(e), 
there is no solution (this region is inaccessible), (ii) For 
h > /^m^rl(e), there is only one stable state 6s > with 
XS < 0. 
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Fig. 33. Magnetization b versus magnetic field h in the micro- 
canonical ensemble for e = — 3 < — 2. 




Fig. 34. Magnetic susceptibility x as a function of the inverse 
temperature 77 in the canonical ensemble for /i — !> 0. 



5.3 The limit h^O 

In this section, we study the magnetic susceptibility in the 
limit ft, — > and investigate its divergence at the critical 
inverse temperature r^c = 1 or critical energy ec = 1 (see 
also [HI] with another presentation). 

Let us first consider the canonical ensemble. For any 
inverse temperature rj < r]c = 1, we know that 6 — 
when /i — > Therefore, we can expand equations 

(Plj) and (j92p in powers of a; ^> at fixed rj. We obtain 
the equivalents 



b{x) 



X 

2' 



h{x) 



V J 2' 



(99) 



from which we deduce the expression of the magnetic sus- 
ceptibility 



X = 



1 — ?/ 



(100) 



This equation is valid for any 77 < ?7c = 1. In particular, 
the magnetic susceptibility diverges at the critical point 

(r; — > ?/,7 = 1^) like 



1 



1 — ?/ 



(101) 



For 77 > = 1 and /i — !■ 0, the magnetic susceptibility 
xiv) is given by equations (|94)) and (|32l) with h = 0. These 
equations can be written 



1 



2ri{x)b'{x) 



T]{x) 



2b{x) 



(102) 



(103) 



For 77 — > 77c = 1, we can expand the previous equations in 
powers of a; — > 0. We obtain the equivalents 



r]{x) - 1 



Xix) 



(104) 
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Fig. 35. Magnetic susceptibility x as a function of the energy 
e in the microcanonical ensemble for /i — ^ 0. The susceptibility 
is negative for e < em- 



so that the magnetic susceptibility diverges at the critical 
point {rj ^ 1]'^ = 1+) like 



1 



X 



2(77-1)- 



(105) 



Note that equations (|10ip and p05p differ by a factor 
2. There is the same difference for the Ising model in a 
magnetic field. The magnetic susceptibility is plotted as a 
function of 77 in Figure 1341 

Let us now consider the microcanonical ensemble. For 
any energy e > Sc = 1, we know that 5 — > when ft — > 
[TT]. Therefore, we can expand equations (^5]) and (ITOj) in 
powers of x — at fixed e. We obtain the equivalents 



6(x)^|, ft(x)^(6-l)|. 



(106) 



from which we deduce the expression of the magnetic sus- 
ceptibility 



1 



X 



e- 1 



(107) 



This equation is valid for any e > Cc = 1- In particular, 
the magnetic susceptibility diverges at the critical point 
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(e — !> e+ = 1+). For e < ec = 1 and h 0, the magnetic 
susceptibility x{e) is given by equations ((98)) . ((32)) and (|33l) 
with /i = 0. These equations can be written 



1 - 2xb{x) 



xb'{x) 



1 + 2x5(a;) ' 



(108) 



(109) 



For e ^ = 1 
in powers of a; — 



' ^ , we can expand the previous equations 
0. We obtain the equivalents 



1 - e(x) 



5x^ 



4 
5x^ ' 



(110) 



so that the magnetic susceptibility diverges at the critical 
point (e — ^ = 1^) like 



1 



X ■ 



2(1 -e) 



(111) 



Note that equations (|107p and piip differ by a factor 2. 
The magnetic susceptibility is plotted as a function of e 
in Figure 1551 

Finally, we compute the behaviour of the magnetiza- 
tion at the critical point when /i — >■ 0. In the canonical 
ensemble, taking 77 = = 1 in equations (|9T|) and (|92|) . 
we obtain h = Ii(x) / Iq(x) and h = x/2 — h{x). Consid- 
ering now the limit /i — ?• 0, corresponding to 2; -> 0, and 
using the asymptotic expansions of Appendix El we ob- 
tain h ^ xjl and h ^ x^/16 yielding 



(2/1) 



1/3 



(112) 



for /i — > at 77 = = 1. This corresponds to the critical 
isotherm. In the microcanonical ensemble, taking e = Ec = 
1 in equations (j95p and (j96p . then considering the limit 
/i — >■ 0, corresponding to .t — > 0, and using the asymptotic 
expansions of Appendix |21 we obtain b ^ x/2 and h ^ 
5x'^/16 yielding 



1/3 



(113) 



for /i— >-Oate = ec = l. The exponent 1/3 is the same in 
the two ensembles. 



heats arc always positive (in canonical and microcanon- 
ical ensembles) but they undergo a discontinuity at the 
critical point (ec,77c) = (1,1) when H = 0. Their mag- 
netic susceptibilities are positive in the canonical ensem- 
ble but they arc negative in the microcanonical ensem- 
ble for emin[h) < e < e„i ~ 0.44025 and they diverge 
at the critical point {ec,ric) = (1,1) when H = 0. We 
have also found metastablc states that arc anti-aligned 
with the magnetic field (BH < 0). These states have 
negative specific heats in the microcanonical ensemble 
for e{rj(.{h)) < e < ec{h) leading to ensembles inequiv- 
alence and zeroth order phase transitions similar to the 
isothermal collapse and the gravothcrmal catastrophe in 
astrophysics. Their magnetic susceptibilities are positive 
in the canonical ensemble but they are negative in the 
microcanonical ensemble for c'^mi^) ^ ^ < when 
< /i < /im ~ 0.2649 and for e^„(/i) < e < edh) when 
hm < h < 1. The magnetic curve B{H) displays there- 
fore hysteretic cycles involving positive or negative mag- 
netic susceptibilities. We have also shown the existence 
of gaps in the magnetization leading to ergodicity break- 
ing. As a result, the HMF model with a magnetic field 
presents a rich phenomenology. It would be interesting to 
compare these theoretical predictions with direct numer- 
ical simulations. Unfortunately, we have shown that the 
"mctastable" states are in fact unstable with respect to 
perturbations that change the phase of the magnetization. 
Since all types of perturbations should be considered in 
principle, this throws doubts on their physical relevance. 
This is a pity because much of the richness of the problem 
disappears. Indeed, if we consider all types of perturba- 
tions, only fully stable states remain and the system does 
not display any phase transition. However, if we consider 
only perturbations that are symmetric with respect to the 
axis determined by the magnetic field H, we should be 
able to observe the phase transitions between aligned (sta- 
ble) and anti-aligned (metastable) states described in this 
paper. 



A Asymptotic expansions 

In this section, we give the asymptotic expansions of b{x), 



ri{x) and e{x) and x(^) for x 
and X x^{h). 

For X- — > 0, we have 



0, 



+ CXD, X 



b{x) 



X 

2 



X 

16 



X 

96 



■o{x^), 



(114) 



6 Conclusion 

In this paper, we have studied the thermodynamics of 
the HMF model with a magnetic field. In the absence 
of magnetic field, the Hamilonian of the HMF model is 
invariant under the translation operation Oi 9i + (j) 
which is equivalent to the C/(l) rotational symmetry. The 
U{1) symmetry is spontaneously broken by the inclusion 
of an external magnetic field H. The fully stable states are 
aligned with the magnetic field {BH > 0). Their specific 



r]{x) 



X X 



2h 4/i2 8h^ 32/i4 



x^ + oix"), (115) 



2h 



eix 



„ , hx 7X , r; , 

1 ~ 2hx 1 1 V oix^) 

8 4 48 ^ ' 



X-cano (•^) 



a; 3x^ 



2h 16h 16h^ 



(116) 



+ oix^), (117) 
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X™(x) = --^-^ + o(x). (118) 



B Absence of critical point with By ^ 



The density profile corresponding to the distribution func- 
For x ^ +CX}, we have tion (HH) is of the form 

6(x) = 1 - - + o(a;-3), (ng) p(e) = Ae^^^- (132) 

Let us write Bx + H = r cos (j) and By = r sin with 
_^ ^ 2 _^ ^ r = [(B:j + H)'^ + Bl] Then, equation becomes 

''^''^ = W+h) + 4(1 + /i)2 + 16(1 + hfx + M « 

(120) ^(^)-2^^"^°^"-''' (1^^) 

where we have used equation ([5]) to determine the am- 
4(1 + /i) 2 — h _3 phtude. Substituting equation (|133p in equations ^ and 

e(x) = -2-4/iH - 2x^ ^ 0, we obtain the self-consistency relations 

27rBx hifir) 
27ri3y /i(/3r) . 

■ , / = ^ , „ , sm 0, (135) 

Xmicroix) ^ - , . + , + 0{X (123) 

z[i -h /ij 0(^1 -t n) X determining the magnetization components and By. 

Assuming By ^ 0, we obtain 



For X — cxD, we have 



Bx + H ^ COS0 ^ 

^ ' -3^ n94^ By sinci " B,/ ^ ^ 



6(a;) = -1 - — + — +o(a;-3), (124) 



2x ' 8x2 



/ ^ ^ I ]^ |_ ^ ^ I „('™-2\ energy or critical point of free energy at fixed mass with 

' 2(1-/1) +4(1-/1)2 16(1- /i)3x"^°^^ By^O. 

(125) 



This relation cannot be satisfied unless H = 0. Therefore, 
there is no critical point of entropy at fixed mass and 



C Eigenvalue equation 

e(x) = -2 + 4/i - 2x2~+'''-^ ^^'^^^ In [H], it is shown that the thermodynamical stability 

problem can be reduced to the study of an eigenvalue 



4(1-/1) 2-/1 , ,_3^ 

equation. The eigenvalue equation associated with the 



1 3 — h 1 , I-, , , maximization problem ( 1131) can be written 

= -2(1^+4(1^^ 



x™(.)--^+^(^+o(^-).(i28) =S^1r+2^^^' (13^) 

For X -> xo{h), we find that ^i*'^ 

/■Sir 

6(x) =6o(/i)-6'(xo(/i))(xo(/i)-x) + ..., (129) V= -^q{e)de, (138) 



^.^(/j-) and q{0) = g(27r) = 0, where q{e) ^ Sp{0')de' is the 

ri{x) ^ ^ 2b'(x (h)](x (h) — x)' (1^^) perturbed integrated density. This corresponds to equa- 

\Xo{ jjy^oK ) X) tions (54) and (55) of [TT] where we have incorporated 

the effect of the magnetic field in the potential {<Ptot = 

2&'fx f/i")") ^ ~'~ ~ —{B + H) cos 9). A critical point of entropy at 

e(x) — eo{h) {xo{h) — x) + ... (131) fixed mass and energy is a maximum iff all the eigenvalues 

^0(^1) A are negative and it is a saddle point if at least one of 
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these eigenvalues is positive. Here, we shall determine the 
point of marginal stability eorresponding to A = 0. We 
follow a method similar to the one developed in [5T] and 



Taking A = in the eigenvalue equation (jl37p . we 
obtain 

T — — oBy COS + oBx sm = 0, 

de \pd9j " MT de 

(139) 



with 



SBy^—[ q{9)cos9d9^~— [ q'{9)smdde, 

(140) 



SB^, = — / q{d) sin 9 d9 = — [ q' {9) cos Ode, 



2tt 



27r 



(141) 



where we have used an integration by parts to obtain the 
seeond equalities. Equation ()139p can be integrated once 
to yield 

^P=^^ ^i^By sin9 + SB., COS9 + C) 
an 1 



2V 



{B + H)pcos9, (142) 



where C is a constant of integration. Another integration 
with the boundary condition q{0) = yields 

q{9)^^-^ f psm9'd9' +^-^ f pcos9'd9' 
^ Jo ^ Jo 

C 2V 
+ - pd9' -j^_{B + H) pcose'd9'.{U3) 

Then, the condition q(2TT) = determines the constant 



C 



-5B, + ^{B + H 



M 



/•2tt 

/ p cos 9 d9, (144) 

"'0 



where we have used the fact that the density profile p{9) is 
symmetric with respect to = to simplify some terms. 
Using an integration by parts, and the expression of 'Ptot, 
the variable V can be written 



V^{H + B) / q'[9)cos9d9. 
Jo 



(145) 



Substituting dq/d9 from equation (jl42p in equation (jl45p 
and solving for V , we obtain 



V 



H + BSB^ J^"" p cos^ 9d9 + C J^"" p cos 9 d9 
~ 1 + j^{B + HY p cos-' 9 d9 



(146) 



Finally, substituting this expression in equation (jl44p and 
solving for C, we get 



C 



SBx p cos 9 d9 



D 



(147) 



where 



2 r^^ 

D = M+—{B + Hfi pcos'9d9 
^ Jo 

;{B + Hf(^J^ pcos9 



MT'' 



2tt \ 2 

d9] . (148) 



Substituting equation (jl42p into equation (|140p . we find 



either that 5By = or that 



1 = 



27TT Ja 



psm^9d9. 



(149) 



Introducing the dimensionless variables defined in Section 
and using the relation 



27r 



X 



(150) 



that can be derived from the density profile ([25|). we find 
that the condition (|149p can be rewritten 



1 = 277- 



(151) 



When = 0, comparing equation (|15ip with equation 
([5^ . we see that this condition is always satisfied. There- 
fore, if we restrict ourselves to perturbations such that 
5Bx = Q and 5 By ^ 0, we conclude that the system is 
always marginally stable j21j . Such perturbations corre- 
spond to a variation of the phase of the magnetization 
vector B. This is just a mere rotation of the equilibrium 
profile. Since the HMF model without magnetic field is 
invariant by rotation, these perturbations do not change 
the entropy. By contrast, when h ^ 0, this U{1) symmetry 



is broken and the condition p5ip can never be fulfilled. 
In fact, the aligned phase is always stable with respect to 
perturbations of the form SB, = and SBy ^ 0, while 
the anti-aligned phase is always unstable by such pertur- 
bations (see Appendix ID)) . There is therefore no marginal 
point corresponding to A = 0. 

Substituting equation (|142p into equation (jl4ip . and 
using the expressions (|146p and (|147p of V and C we find 
either that SB, = or that 



2ttT 



peos 9 d9 



k{B + H) 
^ ttMT^E 

k ^ '•2'^ 



27r ^ 2 

p cos' 9 d9 



2nTDE 



p cos 9 d9 



(152) 



where 



E = 1 



MT' 



(B + H 



•2tt 

i2 / „„„„2 



p cos' 9 d9. (153) 
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Introducing the dimcnsionlcss variables defined in Section 
12.31 and using the relations 



27r 



— / p COS 9d9 ^ b, 



^ M 



2ir 



pcos 



ide = 1 - -, 

X 



(154) 



(155) 



that can be derived from the density profile (j25p . we find 
that the condition (|152p can be rewritten 



1 = 2(7,^x2) 1 



2r;52 



1 + 2a;26'(a 



(156) 



where we have used equations ((32)) and (p9| to simplify 
some terms. We can easily check that r]{x) — = is 
not solution of this equation. Then, after some transfor- 



mations (involving a division by rj - 
equation (|156p is equivalent to 



0), we find that 



[1 - 2(,? - x^)b\x)] 



1 - 



h{x) 



1 

2^ 



= 0. (157) 



We have checked that the last term is strictly positive. 
Therefore, the criterion (|152p is equivalent to 



1 - 2{ri - x^)b' {x) =:0. 



(158) 



For ^ 0, this relation selects the turning point of en- 
ergy ec{h) of the anti-aligned phase (see equation (|H7|)). 
For = 0, the solution of equation (|158p is x = which 
corresponds to the bifurcation point Ec = 1- Therefore, 
the condition of marginal stability in the microcanonical 
ensemble (A = 0) coincides with the turning point of en- 
ergy (if /i 7^ 0) or with the bifurcation point (if h = Q) in 
agreement with the Poineare theorem. 

Let us now consider the canonical ensemble. To that 
purpose, it suffices to take 1^ = in the foregoing expres- 
sions. This yields equation (|149|) and 



2itT 



p cof? 9 d9 



2tvTM 



2w ^2 

p cos 9 d9 



(159) 



This is a particular ease of equations (F.8) and (F.9) in 
|21j . Introducing the dimensionless variables defined in 
Section [231 and using the relations (|151|) . (|155p and (p9|) . 
we obtain 



1 - 2Tjb'{x) = 0. 



(160) 



For h ^ 0, this relation selects the turning point of temper- 
ature r]c{h) of the anti-aligned phase (see equation dM])). 
For h = 0, the solution of equation (jl60p is x = which 
corresponds to the bifurcation point rjc ^ I. Therefore, the 
condition of marginal stability in the canonical ensemble 
(A = 0) coincides with the turning point of temperature 
(if h 7^ 0) or with the bifurcation point (if /i = 0) in 
agreement with the Poineare theorem. 



D Local equivalence of the variational 
problems 

In this Appendix, we study the local equivalence of the 
variational problems and ([5^ in the canonical en- 
semble and of the variational problems (p5|) and ([75)) in 
the microcanonical ensemble. 



D.l Canonical ensemble 

A critical point of (|^5)) is determined by the variational 
principle SF + aTSM = where a is a Lagrange multi- 
plier accounting for the conservation of mass. This leads to 
the distribution ()25p . The magnetization B is obtained by 
substituting equation (|^5)) in equation ([H]) leading to the 
self-consistency relation (|^ . Using the results of |61I66] 
concerning the second variations of free energy, and in- 
troducing the dimensionless variables of Section 12.31 this 
critical point is a (local) minimum of F at fixed mass iff 



6'f = -2ii6b.,r + (Sbyf) + 



{6pf 



d9 > 0, (161) 



for all perturbations Sp that conserve mass: J dpd9 ~ 0. 
The corresponding variations of magnetization are 



M 



M 



6p cos 9 d9, 



Sp sm9 d9. 



(162) 



(163) 



We can always write the perturbations in the form 

Sp = Sp\i+6p^, (164) 

where Sp\\ = [p + Vx cos 9 -\-Vy sin 9)p and 5p±^ = 5p — . 
The second condition ensures that all the perturbations 
are considered. We shall now choose the constants p, 
and Vy such that 



5p\\ d9 = 0, 



1 [ 

6bx = — / Sp\\cos9d9, 
M J ^" 



Sby ^ YT I ^P\\ sin^rf^^ 



M 



This implies that 



Sp± d9 = 0, 



(165) 
(166) 
(167) 

(168) 



6p± cos 9 d9 ^ 0, 



dp^sin9d9 = 0. (169) 
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The conditions (jl65p . (|166|) and (jl67p lead to the relations 
^ + iy^,b 0, (170) 



Sbx = fib + VxJ, 



by ^ Uy{l~l), 



where we have defined 



1= — I pcos^9d0. 



(171) 



(172) 



(173) 



This forms a system of three algebraic equations that de- 
termines the three constants /i, and Vy. Using the equi- 
librium density profile (P5|). and the identity 



2n 

In-l{x) ~ In+l{x) = I,i{x), 

X 



we find that 



1 = 1-- =h\x)+b^ 



(174) 



(175) 



where we have used the equations ([5l]) and (j39[) . Solving 
equations (|170|) . (|17ip and (|172p for /x, z^a; and Vy, and 
using the result (|175p . we find that 



z., - ^5by. (177) 

Therefore, the perturbation 5p\\ takes the form 
1 „. a; 



<5p 



I = -^p{e){b - cose)5b^ + -p{6)sme5by. (178) 



By construction, and are orthogonal for the scalar 
product weighted by 1/p in the sense that 



5p\\5pA 



de = 0. 



(179) 



Indeed, we have 

5p\\5pA 



P 



d6 ~ {p- + Vx cos 9 + Vy sin 0)6 p± dd ~ 0, 



(180) 



where we have used equations (|168p and (jl69p to get the 
last equality. As a result, we obtain 



(Spf 



de 



de 



dO. (181) 



Using equations (|178p and (jl75p . we obtain after simplifi- 
cation 



Therefore, the second order variations of free energy given 
by equation ()16ip can be written 



5^1 



rjb'{x) 



(Sbx)' 



7]b 

J_ f {5PA_? 
r]M 



2 ) {Sbyf 

de. (183) 



Using equations ([5^ and ([5^ . we finally obtain 

5\f-\f"m6bxf^^-^{5byY + ^, I ^-^de. 

(184) 

For the aligned phase {hb > 0), we conclude from this ex- 
pression that p is a (local) minimum of F[p\ at fixed mass 
iff 6 is a (local) minimum of /(5) (the argument is essen- 
tially the same as the one given in IHI])- On the other 
hand; when h 0, the anti-aligned phase {hb < 0) is 
always unstable with respect to (odd) perturbations that 
change the phase of the magnetization, i.e. Sby ^ (it suf- 
fices to take 6bx = and Sp± = 0). Note that, for h = 0, 
this instability disappears due to the U{1) rotational sym- 
metry of the system. Finally, if we impose Sby = 0, i.e. if we 
restrict ourselves to even perturbations, then we see from 
equation (|184p that p is a (local) minimum of F[p] at fixed 
mass iff 6 is a (local) minimum of f(b). Combining these 
results with those of Section[3l we conclude that the states 
(S) are always stable, the states (U) are always unstable 
and the states (M) are metastable for perturbations with 
Sby = while they are unstable for perturbations with 
Sby ^ 0. 

We can also obtain these results by a slightly different 
method inspired by [32]. Considering even perturbations, 
the stability condition ()16ip becomes 



1 



de > 0, 



(185) 



for all perturbations 6p that conserve mass: J SpdO = 0. 
We see that when J SpcosOde = 0, this condition is au- 
tomatically verified. Therefore, we can restrict ourselves 
to perturbations such that / Spcose de ^ 0. On the other 
hand since the functional (5^/ is quadratic in (5p, we can 
impose / Sp cos e de ^ 1 without loss of generality. We now 
look for the perturbation Sp that minimizes (5^/ with the 
constraints J Spde = and J SpcosOdO = 1. Introducing 
Lagrange multipliers p and v^, we find that this pertur- 
bation ia S p = {fi + Vx cos 0)p (it really corresponds to the 
minimum of (5^/ since S^{5'^f) = l/{r]Mp) > 0). Using the 
constraints, we obtain i/x = l/(/— 6^) and p = —b/{I—b'^). 
Finally, substituting this perturbation is equation (jl85p . 
we obtain {6^f)min = ^/ivi^^b^)] ~2. Therefore, inequal- 
ity ([T^ is satisfied iff l/[r]{I - fo^)] - 2 > 0. We can check 
that this condition corresponds to criterion (138) of |22j . 
On the other hand, the equality — &^)] —2 = 

(marginal stability) returns equation (|159p . Finally, us- 
ing equations (|175p and ([CT)) . we find that {5^f)r. 



l/{r,b'{x)) 



\f"{b). Therefore, inequality (fTSSl) is 
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satisfied iff 6 is a minimum of free energy f{b). Consider- 
ing now odd perturbations, the stabihty condition (|16ip 
becomes 



de > 0, 



(186) 



for all perturbations dp. We see that when / Sps'md dd = 
0, this condition is automatically verified. Therefore, us- 
ing the same argument as above, we can restrict our- 
selves to perturbations such that j Spsuvd dO = 1. We 
now look for the perturbation 5p that minimizes (5^/ 
with the constraint J SpainO d0 = 1. Introducing a La- 
grange multiplier fy, we find that this perturbation is 
Sp = VxSuiOp (it really corresponds to the minimum 
of S^f since 6'^{Pf) = l/{r]Mp) > 0). Using the con- 
straint, we obtain Vy = 1/(1 — /). Finally, substituting 
this perturbation is equation (jl86p . we obtain (5"^ f)rnin = 
^ -^)] ~ 2. Therefore, inequality (|186p is satisfied iff 
1/[77(1 — /)] — 2 > 0. We can check that this condition cor- 
responds to criterion (126) of [H]. On the other hand, the 
equality 1/[??(1 — /)] — 2 = (marginal stability) returns 
equation (|149p . Finally, using equations (|175p and (P^ . 
we find that (5^)™*™ = x/{rib) - 2 = 2h/b. Therefore, 
inequality (|186p is satisfied iff bh > 0. 



D.2 Microcanonical ensemble 



A critical point of (|65p is determined by the variational 
principle SS — aSM = 0, where a is a Lagrange multiplier 
accounting for the conservation of mass. This leads to the 
distribution (|25p where the temperature is determined by 
the energy according to equation (j64p . The magnetization 
B is obtained by substituting equation (j25|) in equation 
1^ leading to the self-consistency relation ([55)1 . Using the 
results of |61T66j concerning the second variations of en- 
tropy appropriately generalized to take into account the 
magnetic field (this amounts to replacing by <l>tot), and 
introducing the dimensionless variables of Section [273l this 
critical point is a (local) maximum of S at fixed mass iff 

6^^-!, [ i^de + Tj [{6b^r + {Sbyf] 



M 



-47f{b + hf{Sb,f<0, (187) 



for all perturbations Sp that conserve mass: J SpdO ~ 0. 
We note that the second order variations of entropy p87p 
are related to the second order variations of free energy 
(fTHTll by 



(188) 



Writing the perturbation Sp in the form (jl64p with equa- 
tion ()178p and using equations (jl84p and ((88)) . we obtain 



dB. 



(189) 



From this identity, we arrive at the same type of con- 
clusions as in the canonical ensemble (see the paragraph 
following equation (|184p ). 

We can also proceed as explained at the end of Sec- 
tion ID. II Considering even perturbations, we find that 

{S'^S)r,^ax = - ^Vi^^ f)rmn - W {b + hf = ^s" {b) whcrC we 

have used equation ((88)) to get the last equality. Therefore, 
the system is microcanonically stable with respect to even 
perturbations iff 6 is a maximum of entropy s{b). We can 
check that this condition corresponds to criterion (134) of 
j22j . On the other hand, the equality s"{b) = (marginal 
stability) returns equation (jl58p . Considering odd pertur- 
bations, we find that {6'^s)max = -\'n{5'^ Dmm = -r]h/b. 
Therefore, the system is microcanonically stable with re- 
spect to odd perturbations iff bh > 0. The condition of 
marginal stability {6'^s)max = returns equation (|149p . 



E Partition function and density of states 

In this Appendix, we compute the partition function and 
the density of states of the HMF model with a magnetic 
field and make the connection with the entropy s{b) and 
the free energy f{b) of Sections 0] and |3] (a more general 
discussion is given in Appendix C of [61]). 

In the microcanonical ensemble, the accessible config- 
urations (those having the proper value of energy) are 
equiprobable. Therefore, the probability density of the 
configuration {9i,vi, ...,On,vn) is Pn{Oi,vi, ...,9n,vn) = 
j^6{E — H) where g{E) is the density of states 

g{E) = J S{E - H)d0idvi...deNdvN- (190) 

In other words, g(E)dE gives the number of microstates 
with energy between E and E + dE. The entropy is de- 
fined by S{E) — \ng{E). Integrating over the velocities in 
equation p90p . a classical calculation leads to 

27r^/2 r 

9iE) = TTTlvT / [2(^^ - U)]— d0i...d0Ar, (191) 



r 



where U{9i, ...,9^) is the potential energy (second and 
third terms in the r.h.s. of equation ([1])). Introducing the 
magnetization vector 



k ^ k ^ 



(192) 



i=l 1=1 

The potential energy can be expressed as 



U 



ttB^ 2tt kN , , 



If we impose the constraint X^ili sinOi = (see Remark 
of Section 12. 6p and introduce the dimensionless variables 
of Section [531 the density of states (jl9ip can be rewritten 



gie)^ {e + 2b'+4bh) — n{b)db, (194) 
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where f2(b) denotes the unconditional number of mi- 
crostates corresponding to the macrostate b. The calcu- 
lation of this integral is classical [1]. For N — >• +00, we 
get 



with 



hiX) 



b = 



(195) 



(196) 



Inserting equation (jl95p in equation p94p . we obtain 

g{e) = J e^[^ ln(e+2t^+4h/i)-fcA+ln/o(A)] ^gy) 

Therefore, the density of states (|190p can finally be written 
g{e) = / e^^C-) db, (198) 



where s(6) is given by equation (j78p . In the limit 
+00, we can make the approximation g(e) ~ gW^C'*) 
where 6* is the global maximum of s{b). We finally ob- 
tain limAr_>+tx) j;S{e) = s(6*). 

In the canonical ensemble, the probability 
density of the configuration {9i,vi, ...,0n,vn) is 
PNidi,vi, ...,0n,vn) = zfp)^'^^ where Z{f3) is the 
partition function 

Z{I3)= J e~'^" deidvi...d0NdvN. (199) 

The free energy is defined by F{(3) = — ^lnZ(/3). Inte- 
grating over the velocities in equation (|199l) . we get 



zif3)^(^jj J e-P^de,...de^. 



(200) 



Using equation ()193p . imposing the constraint 
'^i^iSmOi = and introducing the dimensionless 
variables of Section 12.31 and the unconditional number 
of microstates corresponding to the macrostate 6, the 
partition function can be rewritten 



Ziji) = J 6^"^'' +2'''') 12(6) db. (201) 

Inserting equation (jl95p in equation (|20ip . we get 

^(^) =y'e-^[-2b^-46/.+^;,A-^i„/o(A)]^^^ (202) 

Therefore, the partition function can finally be written 

Z{r,) = J e-^f'^^'^ db, (203) 

where /(&) is given by equation ([5^ . In the limit N 
+00, we can make the approximation Z{r]) ^ e^^^^''*') 
where 6* is the global minimum of f{b). We finally obtain 



F Distribution and variance of the 
magnetization 

In this Appendix, we determine the distribution and the 
variance of the magnetization b for the HMF model with 
a magnetic field and show its connection with the free 
energy /(6) and the entropy s{b) of Sections [3] and H] 
The probability density of the magnetization is defined 

by 

P{b) = J S (^Nb - J2 cos o}j S sin 61;^ 

xPN{ei,...,eN)dei...deN, (204) 

where we have imposed the constraint by = X]i!=i sinOi = 
(see Remark of Section 12. 6p and written b for b^ ■ In 
the canonical ensemble, the A'^-body distribution function 
is Pn{0i,vi, ...,0n,vn 
velocities, we obtain 



zf^)^ ■ Integrating over the 



P, 



1 



27r\ 



N/2 



(205) 



where U is the potential energy. In the microcanon- 
ical ensemble, the iV-body distribution function is 
Pn{Oi,vi,...,0n,vn) = ^ H)- Integrating over 

the velocities, a classical calculation gives 



P, 



1 27r^/2 



Recalling that the potential energy can be expressed in 
terms of the magnetization according to equation (jl93|) . 
and introducing the dimensionless variables of Section [^31 
the probability density of the magnetization in the canon- 
ical and microcanonical ensembles is given by 



PcE{b) 



Z{ri) 



(207) 



1 



PMCE{b) = ^(e + 2b'' + Abh)^f2{b), (208) 
5(e) 



where 



f2{b) = J s (^Nb -J2 cos e}j s sin dei...deN, 

(209) 

is the unconditional number of microstates with magne- 
tization b. It is given, for N — >■ +oo, by equation (jl95|) . 
Substituting equation (|195p in equations (|207p and (|208p . 
we obtain 

p^^(6) ^ 1 g-4^[-2fc^-4b/,.+ |6A-2l„/o(A)]^ (210) 

P^jCEib) = 1 e^[i^'^(^+2'''+4^'^)-''^+'"^°(^)]. (211) 
5(e) 
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Fig. 36. Variance of the magnetization as a function of the 
inverse temperature x] in the canonical ensemble. Using the 
resuhs of Appendix^] we find that N {{A\>f) ->■ 1/2 for 77 
(corresponding to a; ^ 0), N{{Abf) ~ 1 / {8{1 + hf r]'^) for -> 
+00 (corresponding to a; ^- +00), N{{Ab)^) ~ 1/(8(1 - hf-q'^) 
for r] +00 (corresponding to x —00) and N{{Ab)^) ~ 
— 1/(2?7) for ri — i- +00 (corresponding to a; — > xo{h)). On the 
other hand, N{{Ab)^) oc ±(77 - r]c{h)y^^^ when 77 77c(/!.) + 
(corresponding to a; ^ x'^°-"°{h)). 
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Fig. 37. Variance of the magnetization as a function of the 
energy e in the microcanonical ensemble. Using the results of 
AppendixEl we find that N{{Abf) 1/2 for e +00 (cor- 
responding to a; -> 0), N{{Abf) ~ (e - e™,„(/i))V(64(l + hf) 
for e — ^ ^min{h) (corresponding to a; —J- -l-cxo), N{{Ab)'^) ~ 
(e - €'^,Jh)f/{6A{l - hf) for e eLn^i) (correspond- 
ing to a: ^ -00) and N{{Abf) ~ -(e - eo(^))/2 for e ^ 
eoC*) (corresponding to a; — xo{h)). On the other hand, 
A'^((Zi6)'^) oc ±{ec{h) — e)"^''^ when e — ec(h)^ (corresponding 
to a; ^ a;r""°(/i)). 



The distribution of the magnetization can therefore be 
rewritten 



Pcsib) 



1 



ziv) 



e 2 



PhlCEib) = —r-re 
5(e) 



JVs(&) 



(212) 



(213) 



where f(b) and s(6) are the free energy and the entropy 
defined in Sections [3] and ID 

For N +00, P{h) is strongly peaked around the min- 
imum of f{h) or around the maximum of s(6). Therefore, 
we can make a Gaussian approximation 



P{b) 



1 



where the variance is given by 
N{{Ahf)cE 



-JAiiL- 



,?/"(fo)' 



N{{Ahf)MCE 



1 



s"(&) 



Using equations ([5^ and ([53]), we obtain 

1 



iV((zi6)2)MC£ = 



2?7' 



(214) 



(215) 



(216) 



(217) 



(218) 



These quantities are plotted in Figs.[Mland[571 In the anti- 
aUgned phase, the variance diverges at the critical inverse 
temperature rjc{h) in the canonical ensemble (see equation 
(|6ip ) and at the critical energy ec{h) in the microcanonical 
ensemble (see equation (|57)) ). Of course, the variance of 
the unstable states is negative. These results generalize 
those obtained in [61] for h = Q. 

Comparing equations and (|217p . we find that 

Xcano = 2?7A^((Z\&)^)c£; which corresponds to the gen- 
eral identity x = ^Pii'^B)^) valid in the canonical en- 
semble. On the other hand, comparing equations (j98p and 
(|218p and evaluating bdb/de from equations ([5^ and ([55]) . 

we find that \rnicro 

~ 2'qN{{AbY)MCE + ^bdb/de which 
corresponds to the general identity x = ^{Pii^B)'^) + 
B{dB I dE)) valid in the microcanonical ensemble for N 
~\-oo (62j . Note that our approach provides an alternative 
derivation of this relation. 



G Distribution and variance of the energy in 
the canonical ensemble 

The distribution of energy in the canonical ensemble is 
defined by 

p{E)^ J s{E-H{ei,vi,...,eN,VN)) 

y.PN{6i-,vi, On, vn) d9idvi...d6NdvN ■, (219) 



where Pm{Oi,vi, ...,9n,vn) 



is the canonical 



A^-body distribution function. Inserting this distribution 
in equation (|219p and introducing the density of states 
dUni), we obtain 



P{E) 



1 



Z{I3) 



9{E)e 



-PE 



(220) 
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Introducing the microcanonical entropy S{E) = In g{E) 
and the diniensionless variables of Section 12.31 we get 



Pie) 



„JV(s(e)-i,,e) 



ziv) 



(221) 



Using the results of Section 12.31 the microcanonical en- 
tropy s(e) is given, in the limit N — > +oo. by 



2{bix) + h) 



+ ln/o(.T) -a;6(a;), (222) 



e = ^^^^"^^ - 2b{xf - m{x), (223) 



6 = 



h{x)_ 



(224) 



The most probable energy e at inverse temperature rj is 
solution of s'(e) = ^77. Using equations (|222D - (|224D . we 
find that it satisfies the relation 

X 



V 



2{b{x) + h)' 



(225) 



This returns equation p2p . Expanding equation (|22ip 
around this most probable energy, we obtain 



P{e) 



1 '•^'K. 

.g 2{(/le)^> ^ 



where the variance of energy is given by 

1 



.s"(e)- 



(226) 



(227) 



Using equations (j222p - (j224p . we can compute s"(e) at 
fixed rj. If we evaluate this expression at the most prob- 
able energy satisfying equation (|225p and substitute the 
result in equation (|227p . we find that 



NiiAef) 



2 l-2&'(x)(?7-a;2) 
^ 1 - 2T]b'{x) ■ 



(228) 



Comparing equation (j228p with equation (|4T|). we obtain 



NiiAef) = -c, 



(229) 



which returns the general relation C = /3^((Z\£')^), where 
C = —(i'^dE/dp, valid in the canonical ensemble. 



H Specific heat in canonical and 
microcanonical ensembles 

We consider a general Hamitonian 

N 2 

if = ^m^-fC/(ri,...,rw) = A- + [/, (230) 

1=1 



in d dimensions. In the canonical ensemble, the partition 
function is 



Z{I3)^ J e'^" dri...drNdvi...dvN, 



(231) 



where /3 = 1/T is the inverse temperature (we take the 
Boltzmann constant equal to unity). The average energy 
is given by 

(E) ^ ^ J He-'^" dYi...drMdwi...dwN 



1 dZ 



(232) 



where F 
have 



dp 813' 
-(l//3)lnZ is the free energy. Similarly, we 



= ^ j H^e-f^" dri...dTNdwi...dwN = 



1 d'^Z 



z 9/32 



(233) 



The specific heat is defined by 
C_d{E)_ ^,d{E) 



& 



P l^iPF)- (234) 



Substituting equation (|232p in equation (|234p we obtain 



1 d^Z 1 fdZ 



Z 9/32 ^2 I Qp 



(235) 



Using equations (|232p and (|233p , and introducing the vari- 
ance of the energy {{AEY) = {E'^) — (E)^ , we obtain the 
well-known formula 



C = l3^{{AEf), 



(236) 



which shows in particular that the specific heat is positive 
in the canonical ensemble. Note that the average value 
of the kinetic energy and its fluctuations in the canonical 
ensemble are given by 



dN 



(237) 



These relations can easily be obtained from the previous 
ones by using the fact that the variables of velocity and 
position factorize (this is equivalent to taking J7 = in 
the foregoing formulae). 

In the microcanonical ensemble, the density of state is 

g{E)^ j 5{E - H)dvi...dYNdwi...dwN. (238) 
Integration over the velocities yields 



9{E) 



dN 

TT 2 



(E -U) — -^ dri...dr 



N- 



(239) 
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Introducing the entropy S{E) ~ lng(_E), the microcanon- 
ical inverse temperature (3{E) = 1/T{E) is given by 



/5 = 



dS_ _ 1 dg 
dE ^ g{E) dE 

Using equation (|239p . we obtain 



(240) 



1^ J{E- U)-^Pn dT^...dvN, (241) 



where P;v(ri, rjv) is the microcanonical distribution. 
Noting that E — U ^ K , this relation can be rewritten 



dN 



(242) 



and it provides an exact relationship between the inverse 
microcanonical temperature and the average inverse ki- 
netic energy. For +oo, we have 



dN 
2^' 



(243) 



which can be compared with equation (|237[) . The specific 
heat is defined by 



C 



dE 
df 



dE 



S"{E) 



Using equation (j240p . we obtain 



1 

C 



1 



1 0^9 
gdE^ 



1 

7 



dg 
dE 



(244) 



(245) 



Using equations (|239p and (|240p . the foregoing relation 
can be rewritten 



1 

C 



1 



dN 



1 



dN 



1 



(246) 



Introducing the variance of the inverse kinetic energy 

2\ / 1 \ / 1 \ 2 

(247) 



4 



1 



the specific heat per particle c = C/N, and using equa- 
tion (|242p . we obtain after some rearrangements the exact 
relation 



N 



dN 
2 



dN 
2 



1 /4 



For N 



-oo, it reduces to 



N((Aj,] 



(248) 



(249) 



This returns the relationship obtained by Lebowitz et al 
[70] . It clearly shows that the specific heat is not neces- 
sarily positive in the microcanonical ensemble. 



I Magnetic susceptibility in canonical and 
microcanonical ensembles 



Using equation (|192p . the Hamitonian ((T)) can be written 

(250) 



n = Ho- ^BH, 
k 



where Hq is the Hamiltonian of the HMF model without 
magnetic field. In the canonical ensemble, the partition 
function is 



z{i3,H) = J e-mo~^BH) dei...deNdvi...dvN. 



(251) 



The average magnetization is given by 

(B) = ^ J Be-'^(^«-TS^) d9i...deNdvi...dv 



N 



k 1 dZ 



k dlnZ 



27rP Z dH 27r/3 dH 



k dF 
2^aff ' 



(252) 



where F 
have 



-(l//3)lnZ is the free energy. Similarly, we 



(£2) = ly s2g-/J(«o-^B//) dei...d9Ndv,...dvN 



( kVld^Z 



The magnetic susceptibility is defined by 
d{B) _ k d^F 



X 



dH 27r 9i72 " 
Substituting equation (|252p in equation (|254p we get 



(254) 



C 



27r/3 



1 d'^Z 1 f dZ 



ZdH^ Z2 \dH 



(255) 



Using equations (j252|) and (j253|) , and introducing the vari- 
ance of the magnetization ((Z\_B)^) = (S^) — (S)^, we 
obtain the exact relation 



(256) 



which shows in particular that the magnetic susceptibility 
is positive in the canonical ensemble. 

In the microcanonical ensemble, the density of states 

is 



g{E,H)=A / [E-Uo + ^BH 



dei...de 



N, 



(257) 
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where Uq is the potential energy of the HMF model with- 
out magnetic field and A = (27r)~/P(Y)- The average 
magnetization is given by 



(258) 



Taking its derivative with respect to H and E, we obtain 
d{B) (N \2n/B^\ Idg 



d{B) 
dE 



N 



-(I 



ldg_ 
~gdE 



(B), (260) 



where K = E — Uq + ^BH is the kinetic energy. On the 
other hand, taking the derivative of the density of states 
(|257p with respect to H and E, we get 



gdH V 2 ) k \K 



gdE \2 )\K 



The magnetic susceptibility is defined by 

d{B) 



X 



dH 



(261) 



(262) 



(263) 



Combining the preceding equations, we obtain the exact 
relation 



kx 
27r 



N 



- 1 



-{B)^. (264) 



dE 



It generalizes the identity 



(265) 



valid in the microcanonical ensemble for TV +oo (see 
Appendix [F| . These relations clearly show that the mag- 
netic susceptibility is not necessarily positive in the mi- 
crocanonical ensemble. 
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